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Foreword  and  Credits 


This  report  was  prepared  at  the  request  of  the  Hydrologic  and  River  Engineering  Section  of  the 
Portland  District,  U.S.  Army  Corps  of  Engineers.  The  HEC  was  asked  to  assist  in  the  assessment  of 
hypothetical  mudflow  events  that  might  occur  if  the  debris  blockage  presently  containing  Castle  Lake  near 
Mount  St.  Helens  WA  were  to  fail.  The  Portland  District  Corps  of  Engineers,  at  the  request  of  the  U.S. 
Forest  Service  (managers  of  the  Mount  St.  Helens  National  Volcanic  Monument),  undertook  their  initial 
studies  in  1988  to  analyze  the  existing  conditions  of  the  blockage,  determine  the  degree  of  risk  for 
downstream  flooding  posed  by  Castle  Lake,  evaluate  alternatives  to  reduce  that  risk,  and  recommend  a 
solution  for  reducing  the  risk.  The  primary  purpose  of  the  present  study  was  to  estimate  the  potential  for 
flooding  downstream  from  the  Corps’  Sediment  Retention  Structure  (SRS)  for  various  hypothetical  lake 
breaching  scenarios.  It  is  not  the  intent  of  this  investigation  to  evaluate  any  aspect  of  the  risk  of  failure. 

It  merely  quantifies  the  downstream  flood  potential  for  various  hypothetical  breaching  scenarios.  Results 
from  this  investigation  arc  to  be  used  by  USFS  and  USACE  managers  to  decide  what  alternatives  may  be 
effective  in  reducing  the  flooding  potential  in  communities  downstream  from  the  SRS. 

Information  and  data  used  during  this  investigation  and  presented  in  this  report  were  obtained  from  the 
Portland  District  Corps  of  Engineers,  the  U.S.  Forest  Service,  the  U.S.  Geological  Survey  and  the 
Washington  State  Department  of  Ecology.  An  annotated  bibliography  of  reference  materials  examined 
during  this  investigation  is  provided  in  Appendix  A. 

Robert  C.  MacArthur  (Principal  investigator),  Gary  Brunner  and  Doug  Hamilton  conducted  the 
investigation  and  wrote  this  report.  Messrs.  MacArthur  and  Brunner  work  at  the  Hydrologic  Engineering 
Center  in  Davis,  CA.  Doug  Hamilton  is  a  principal  engineer  with  RIVERTECH,  Inc.  in  Laguna  Hills,  CA 
and  provided  technical  assistance  to  the  Hydrologic  Engineering  Center  during  this  investigation.  Ron 
Mason  was  the  Project  Manager  for  the  Portland  District  Corps  of  Engineers  and  Mr.  John  Steward  was 
the  Project  Manager  for  the  U.S.  Forest  Service.  Vernon  Bonner  was  the  Chief  of  the  Training  Division 
during  the  study  and  Mr.  Darryl  Davis  was  the  Director  of  the  Hydrologic  Engineering  Center  during  the 
investigation. 

The  investigation  reported  herein  was  conducted  by  the  Hydrologic  Engineering  Center  in  Davis, 
California  at  the  request  of  the  Portland  District,  Corps  of  Engineers  and  the  Gifford  Pinchot  National 
Forest,  U.S.  Forest  Service.  Funding  for  this  study  was  provided  by  the  U.S.  Forest  Service. 


Executive  Summaty 


The  eruption  of  Mount  St.  Helens  on  May  18,  1980  produced  a  debris  avalanche  that  moved  down  tin: 
North  Fork  Toutle  River  damming  several  tributary  streams.  The  blockage  at  the  confluence  of  Soutn 
Fork  Castle  Creek  and  Castle  Creek  produced  a  natural  debris  dam  approximately  190  feet  high.  Snow 
melt  and  runoff  waters  captured  behind  the  blockage  quickly  formed  a  lake.  Even  though  the  U.S.  Armv 
Corps  of  Engineers  installed  an  emergency  spillway  in  October  1981  to  prevent  overtopping,  there  is 
concern  that  the  debris  dam  may  be  unstable  during  the  combination  of  a  severe  hydrologic  event  and  an 
earth  quake  of  magnitude  6.8  or  greater. 

The  purpose  of  this  study  is  to  evaluate  the  hydraulic  characteristics  of  mudflow  events  resulting  from 
the  hypothetical  failure  of  Castle  Lake  and  to  examine  the  ability  of  the  Sediment  Retention  Structure 
(SRS)  to  capture  and  pass  such  events  through  its  emergency  spillway  for  various  initial  conditions  at 
Castle  Lake  and  in  the  SRS.  More  specifically,  the  study  is  to:  (1)  determine  if  flows  will  exceed  the 
present  spillway  capacity  of  the  SRS,  (2)  determine  if  the  SRS  will  be  overtopped  during  various  breaching 
scenarios,  (3)  estimate  how  the  peak  discharge  in  communities  downstream  from  the  SRS  will  be  affected 
by  the  presence  of  the  SRS,  (4)  evaluate  the  routing  effects  on  the  resulting  mudflow  hydrographs  due  to 
lowering  the  initial  Castle  Lake  levels  at  the  time  of  breaching,  and  (5)  evaluate  the  performance  of  the 
SRS  during  these  various  events  when  the  SRS  is  empty  of  water  and  sediment  (existing  conditions),  or 
full  of  sediment  deposits  up  to  the  spillway  crest. 

It  is  not  the  intent  of  this  investigation  to  evaluate  any  aspect  of  the  risk  of  failure.  It  merely  quantifies 
the  downstream  flood  potential  for  various  hypothetical  breaching  scenarios.  Results  from  this 
investigation  are  to  be  used  by  the  USFS  and  USACE  managers  to  decide  what  alternatives  may  be 
effective  in  reducing  the  flood  potential  in  those  communities  downstream  from  the  SRS. 

The  Hydrologic  Engineering  Center  (HEC)  reviewed  available  scientific  and  engineering  literature 
pertaining  to  past  debris  blockage  failures,  the  characteristics  of  flow  bulking,  and  methods  for  simulating 
the  breaching  of  the  dam  and  routing  of  the  resulting  dambreak  hydrograph  downstream  to  the  SRS  and 
beyond  to  the  Columbia  River.  HEC  examined  more  than  half  a  dozen  different  breaching  models  and 
decided  to  use  the  National  Weather  Service’s  BREACH  model  (Fread,  1989)  because  it  is  most  complete 
model  available.  It  is  also  physically  based  with  respect  to  the  development  of  different  modes  of  dam 
failure.  BREACH  uses  soil  properties,  sediment  transport  functions  and  hydraulic  computations  to  predict 
the  breaching  characteristics  and  discharge  hydrograph  emanating  from  a  breaching  earthen  dam  or  debris 
blockage.  The  BREACH  model  was  used  to  develop  breach  outflow  hydrographs  for  several  types  of 
breaching  scenarios  and  various  lake  levels. 

Hydrographs  were  developed  for  three  different  initial  water  surface  elevations  in  Castle  Lake:  (1)  2,580 
feet  above  NGVD,  (2)  the  lake  lowered  30  feel  to  2,550  NGVD  and  (3)  the  lake  lowered  60  feet  to  2,520 
NGVD.  The  National  Weather  Service’s  DAMBRK  model  (Fread,  1989)  was  used  to  route  the 
damhrcach  hydrographs  down  the  valley,  through  the  Sediment  Retention  Structure  (SRS),  all  the  way 
downstream  to  the  Columbia  River.  Energy  based  procedures  were  developed  to  simulate  the  bulking  of 
the  flows  via  a  series  of  lateral  inflow  hydrographs.  The  hydrographs  were  shaped  and  positioned  along 
the  routing  reach  to  provide  the  appropriate  timing  and  volume  of  the  lateral  inflow  according  to  the 
characteristics  of  the  primary  flood  wave  in  the  channel. 

Dambrcach  hydrographs  and  associated  lateral  bulking  hydrographs  were  developed  for  three  initial  lake- 
levels  in  Castle  Lake  and  for  two  breaching  scenarios:  (1)  a  piping  failure  due  to  "heave"  as  per  the 
USGS’s  report  by  Lacnen  and  Orzol  (1987)  with  a  modification  of  the  flow  bulking  as  per  HEC  method-, 
and  (2)  a  piping  failure  positioned  over  the  historical  South  Fork  Castle  Creek  outlet  channel  (referred  to 
as  the  HEC  or  Corps  Recommended,  Breaching  Scenario).  Downstream  bulking  of  the  flows  depends  on 
the  initial  volume  and  duration  of  the  outflow  hydrograph,  on  the  breaching  mechanisms  and  on  the  valle\ 
soil  properties  and  water  content.  Soil  samples  were  collected  from  the  downstream  valley  debris  deposits 
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by  the  USGS  (Meyer  and  Dodge,  1^88)  and  the  Corps  of  Engineers  (l ’SACK  O' l  and  <!.»!.■ 

USACE,  1990).  A  range  of  measured  values  for  the  key  parameters  used  to  determine  bulking  .»<•« 
mudflow  characteristics,  such  as  porosity,  percent  saturation,  and  expected  sediment  concentrations,  was 
developed.  A  Monte  Carlo  weighting  technique  was  utilized  to  determine  the  most  probable  combination 
of  these  parameters.  From  the  results  of  the  Monte  Carlo  simulations,  high,  medium,  and  low  hulking 
factors  were  selected  as  a  range  of  probable  values  for  (he  sensitivity  analysis  that  was  conducted  bv  MFC' 

Final  breaching  and  routing  simulations  were  conducted  based  on  what  is  referred  to  throughout  the 
remainder  of  the  report  as  "the  HE’C  or  Corps  Recommended,  Breaching  and  Bulking  Scenarios."  They 
represent  the  breaching  and  bulking  characteristics  recommended  and  agreed  upon  by  the  Corps  of 
Engineers  and  the  U.S.  Forest  Service.  'Hie  estimated  peak  discharge  from  a  hypothetical  failure  of  the 
Castle  Lake  blockage  (using  the  Corps’  breaching  and  hulking  scenario)  exceeds  the  peak  discharges 
predicted  from  potential  energy  versus  peak  discharge  relationships  developed  from  historical  dam  failures 
by  more  than  2.3  times.  It  exceeds  the  predicted  peak  discharge  envelope  curve  from  historical  dam 
failures  by  3.6  times.  Therefore,  the  Corps  recommended  breaching  and  bulking  scenario  produces  a  larger 
discharge  (i.c.,  flows  that  are  larger  than  the  historical  envelope  curves)  because  the  breaching  scenario 
depicted  represents  a  simultaneous  occurrence  of  a  severe  hydrologic  event  and  an  earthquake  of 
magnitude  6.8  or  greater.  The  resulting  dam  failure  will  produce  a  greater  Hood  hydrograph  than  a 
traditional  overtopping  type  of  failure  analysis. 

Routing  results  indicate  that  the  SRS  significantly  reduces  the  peak  discharge  from  the  hypothetical 
failure  of  Castle  Lake.  The  initial  elevation  of  Castle  lake  prior  to  failure  also  affects  the  magnitude  ol 
the  resulting  dambreak  discharge.  For  the  failure  and  bulking  scenario  recommended  by  the  Corps  ol 
Engineers,  the  SRS  reduces  the  peak  discharge  into  the  North  Fork  Toutle  River  by  85  percent  (from 
695,000  to  105,200  cfs)  for  full  lake  conditions.  If  Castle  Lake  is  lowered  30  feet  prior  to  its  failure,  the 
SRS  reduces  the  peak  flow  by  82  percent  (from  352, 500  to  62,000  cfs).  If  Castle  Lake  is  lowered  by  60 
feet,  the  SRS  reduces  the  peak  flow  by  95  percent  (from  131 ,8(K)  to  6, (XX)  cfs).  The  amount  of  storage  the 
SRS  can  provide  depends  on  how  full  of  sediment  it  is  when  a  flood  event  occurs.  Under  the  present 
"existing  conditions"  in  the  SRS  and  the  Cowlitz  River,  all  of  the  Corps  recommended  flooding  scenarios 
would  be  fully  contained  within  the  channel  at  Castle  Rock  and  Kelso  -  Longview.  The  resulting  Hows 
would  be  similar  to  a  1%  chance  flood  event  in  the  Cowlitz  River.  If  the  SRS  were  full  of  sediment,  all  of 
the  Corps  recommended  flooding  scenarios  would  be  fully  contained  within  the  channel  at  Kelso  - 
Longview,  but  not  at  Castle  Rock  for  cither  lake  full  or  lake  lowered  30  feet  conditions.  With  the  SRS 
initially  full  of  sediment,  it  is  estimated  that  the  channel  eapacity  near  Castle  Rock  will  be  exceeded  unless 
the  initial  lake  elevation  is  lowered  by  60  feet.  None  of  the  hypothesized  breaching  and  bulking  scenarios 
will  exceed  or  overtop  the  SRS  for  cither  "existing  conditions"  or  "full  conditions." 

The  SRS  protects  communities  and  those  river  sections  downstream  from  it  but  docs  not  affect  the 
areas  upstream  from  the  SRS.  Therefore,  there  is  no  protection  above  the  SRS  from  a  hypothetical  failure 
of  Castle  Lake.  Because  of  this,  and  the  possibility  that  a  failure  of  Castle  Lake  would  greatly  reduce  the 
sediment  storage  capacity  and  effective  life  of  the  SRS,  additional  field  investigations  pertaining  to  the 
geotechnical  stability  of  the  Castle  Kike  blockage  and  regular  monitoring  are  recommended. 

Tlie  Castle  Lake  blockage  was  10  years  old  in  May,  1990  and  appears  to  be  stable  under  its  past  and 
present  conditions.  The  Portland  District  Corps  of  Engineers  installed  an  emergency  spillway  in  October 
of  1981  to  stabilize  the  lake  elevation  at  2577  feet  above  NGVD.  Groundwater  levels  in  the  blockage  and 
seeps  along  the  downstream  face  of  the  blockage  have  been  monitored  since  the  eruption.  According  to 
the  Geotechnical  Branch  of  the  Corps’  Portland  District,  there  is  no  field  evidence  of  unstable  conditions 
in  the  blockage  materials  since  the  installation  of  the  spillway.  Even  though  debris  blockage  dams  form  in 
a  wide  variety  of  physiographic  settings,  most  debris  blockage  dams  arc  very  short  lived.  Investigations  of 
historical  landslide  and  debris  blockage  dam  failures  indicates  (hat  approximately  22  percent  of  the 
landslide  dams  failed  in  less  than  1  day  after  formation  and  that  half  failed  within  a  period  of  10  days. 

Less  than  10  percent  of  the  natural  debris  blockage  dams  last  more  than  1  year.  It  was  also  observed  that 
more  than  50  percent  of  the  documented  debris  and  landslide  dams  failed  due  to  overtopping.  In  the  case 
of  Castle  Lake,  the  possibility  for  overtopping  has  been  eliminated. 
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Numerical  Simulation  of  Mudflows  from  Hypothetic! 
Failures  the  Castle  Lake  Debris  Blockage 
Near  Mount  St.  Helens,  \A,X 


1.  Study  Purpose 

The  May  18,  1980  eruption  of  Mount  St.  Helens,  WA,  produced  a  debris  avalanche  that  flowed  down 
the  North  Fork  Toutle  River  damming  several  tributary  streams.  The  blockage  at  the  confluence  of  South 
Fork  Castle  Creek  and  Castle  Creek  produced  a  natural  debris  dam  approximately  190  feet  high.  Figure  1 
shows  the  general  study  area  near  Mount  St.  Helens  and  the  location  of  Castle  Lake.  Snow  melt  and 
runoff  waters  captured  behind  the  blockage  quickly  formed  a  lake.  To  prevent  overtopping  and  a 
potentially  catastrophic  failure  of  the  blockage  retaining  Castle  Lake,  the  U.S.  Army  Corps  of  Engineers 
(USACE)  constructed  an  emergency  spillway  in  October  1981  at  the  eastern  end  of  the  blockage  to 
stabilize  the  lake  at  elevation  2,577  feet  NGVD.  Studies  by  the  U.S.  Geological  Survey  (USGS)  indicated 
that  "the  blockage  is  potentially  unstable  against  failure  from  piping  due  to  heave  and  internal  erosion 
when  groundwater  levels  are  seasonally  high"  and  that  an  earthquake  of  6.8  or  greater  might  initiate  uch  a 
failure  (Lacnen  and  Orzol,  19*7).  If  the  Castle  Lake  blockage  were  to  fail  rapidly  by  the  mechanism 
suggested  by  the  USGS,  approximately  18,500  acre-feet  (AF)  of  stored  water  in  the  lake  could  create  a 
mudflow  flood  event  in  the  North  Fork  Toutle  River.  The  USGS  (Laencn  and  Orzol,  1987)  estimates  that 
an  event  of  this  nature  could  result  in  a  peak  discharge  of  2,100,000  cfs  at  the  Corps’  N-l  debris  retention 
dam  ten  miles  downstream  from  Castle  Lake  (see  Figure  l)  and  possibly  lead  to  downstream  flooding. 

In  the  wake  of  the  Mount  St.  Helens  eruption,  the  Corps  developed  a  long-term  flood  control  and 
navigation  maintenance  plan.  A  major  component  of  that  plan  is  the  $56.5  million  Sediment  Retention 
Structure  (SRS)  designed  to  trap  the  huge  amounts  of  sediment  expected  to  continue  to  move  dow  n  the 
North  Fork  Toutle  River.  The  SRS  was  designed  to  capture  runoff-induced  sediment  from  the  blast  zone, 
thus  preventing  sediment  deposition  in  the  Cowlitz  and  Columbia  Rivers.  Without  the  SRS,  sediment 
materials  could  continue  to  accumulate  in  the  rivers  below  the  SRS  thus  reducing  their  flood  routing  and 
navigation  capacities.  There  was  additional  concern  that  failure  of  the  Castle  Lake  blockage  resulting  in 
the  possible  occurrence  of  a  mudflow  event  could  jeopardize  the  safety  and  performance  of  the  SRS  or 
perhaps  lead  to  flooding  in  communities  downstream  from  the  SRS. 

The  purpose  of  this  study  is  to  evaluate  the  hydraulic  characteristics  of  mudflow  events  resulting  from 
the  hypothetical  failure  of  Castle  Lake  and  to  examine  the  ability  of  the  SRS  to  capture  and  pass  such 
events  through  its  spillway  for  various  initial  conditions  at  Castle  Lake  and  in  the  SRS.  More  specifically, 
the  study  is  to:  (1)  determine  if  flows  will  exceed  the  present  spillway  capacity  of  the  SRS,  (2)  determine 
if  the  SRS  will  be  overtopped  during  various  breaching  scenarios,  (3)  estimate  how  the  peak  discharge  in 
communities  downstream  from  *hc  SRS  will  be  affected  by  the  presence  of  the  SRS,  (4)  evaluate  the 
ro-.iting  effects  on  the  resulting  mudflow  hydrographs  due  to  lowering  the  initial  Castle  Lake  levels  at  the 
time  of  breaching,  and  (5)  evaluate  the  performance  of  the  SRS  during  these  various  events  when  the  SRS 
is  empty  of  water  and  sediment  (existing  conditions),  or  full  of  sediment  deposits  up  to  the  spillway  crest. 
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2.  Approach 

The  Hydrologic  Engineering  Center  (HEC)  conducted  this  investigation  in  two  phases.  The  first  phase, 
a  reconnaissance  level  investigation,  included  a  field  inspection  of  the  Mount  St.  Helens  National  Volcanic 
Monument,  the  Castle  Lake  debris  blockage  area  and  valley  sections  downstream  from  the  blockage  all  the 
way  to  the  SRS.  HEC  staff  attended  two  days  of  meetings  with  project  personnel  from  the  Portland 
District  USACE  and  the  U.S.  Forest  Service  to  discuss  the  background  of  the  problem,  concerns  they  and 
other  Federal  and  State  agencies  had  for  the  safety  of  the  blockage,  and  to  outline  technical  procedures  for 
conducting  the  analytical  investigation.  The  Phase  1  investigation  also  included  a  thorough  literature 
investigation  and  a  "reconnaissance-lever  (preliminary)  mudflow  routing  investigation  (see  Appendices  A 
and  B  respectively).  Appendix  B  contains  a  summary  paper  ihat  presents  the  results  from  the  Phase  1 
Reconnaissance-Level  investigation.  Results  from  the  Phase  1  studies  were  presented  to  project  managers 
from  the  1 JSGS,  USFS,  State  of  Washington  Department  of  Ecology  and  Dam  Safety  (SWDE)  and  the 
Portland  District  IJSACE.  The  results  were  used  to  formulate  an  agreed-upon  analytical  approach  and 
ranges  of  breaching  and  mudflow  bulking  parameters  to  be  used  during  the  Phase  2  studies.  The 
remainder  of  (bis  report  concentrates  on  the  procedures  and  results  from  the  Phase  2  investigation. 

The  Phase  2  investigation  included  the  development  of  energy  based  procedures  for  bulking  and 
debulkmg  the  dam  break  flows  from  hypothetical  breaching  of  the  Castle  Lake  blockage.  The  National 
Weather  Service’s  BREACH  model  (Fread,  1989)  was  used  to  develop  breach  outflow  hydrographs  for 
several  types  of  breaching  scenarios,  and  various  lake  levels.  BREACH  is  a  physically  based  model  that 
uses  soil  properties,  sediment  transport  functions,  and  hydraulic  computations  to  predict  the  breach 
characteristics  and  the  discharge  hydrograph  emanating  from  a  breaching  earthen  dam  or  debris  blockage. 
The  critical  breaching  time  (defined  as  the  time  from  the  beginning  of  a  major  rise  in  the  outflow 
hydrograph  until  the  time  of  the  peak  flow  out  of  the  breach)  for  piping  and  heave  failures  was  determined 
to  be  approximately  15  minutes.  Hydrographs  were  developed  for  three  different  initial  water  surface 
elevations  in  Castle  Lake:  (1)  2,580  feel  above  NGVD,  (2)  the  lake  lowered  30  feet  to  2.550  NGVD  and 
(3)  the  lake  lowered  60  feel  to  2,520  NGVD.  The  National  Weather  Service’s  DAMBRK  model  (Fread 
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19X9)  whs  used  to  route  the  dambrcach  hydrographs  down  valley,  through  the  Sediment  Rc.er.in  > 
Structure  (SRS),  and  continuing  down  to  the  Columbia  river.  Energy  ba-cd  procedures  were  d.  elop*  >J  to 
simulate  the  bulking  up  of  the  flows  via  a  series  of  lateral  inflow  hydtographs.  The  hvdu.t were 
shaped  and  positioned  along  the  routing  reach  so  as  to  provide  the  appropriate  timing  and  vo.ume  of  he 
lateral  inflow  according  to  the  magnitude  of  the  primary  flood  wave  in  the  channel.  Breakout  hvdrogMphs 
and  associated  lateral  bulking  hydrographs  were  developed  for  three  initial  lake  levels  in  C'asfle  I  jkc  and 
for  two  different  breaching  scenarios:  (1)  a  piping  failure  due  to  heave  as  per  the  USGS’s  report  by 
Laencn  and  Orzol  (1987),  and  (2)  a  piping  failure  positioned  over  the  historical  South  Fork  Castle  Creek 
outlet  channel  (referred  to  as  the  HEC  Breaching  Scenario).  Downstream  bulking  of  the  flows  depends  on 
the  initial  volume  and  duration  of  the  outflow  hydrograph,  on  the  breaching  mechanisms  and  on  the  valley 
soil  properties  and  water  content.  Soil  samples  were  collected  from  the  downstream  valley  debris  deposits 
by  the  USGS  (Meyer  and  Dodge,  19X8)  and  the  Corps  of  Engineers  (USAGE,  1984  and  unpublished  data, 
USAGE,  1990).  A  range  of  measured  values  for  the  key  parameters  used  to  determine  bulking  and 
mudflow  characteristics,  such  as  porosity,  percent  saturation,  and  expected  sediment  concentrations,  were 
developed.  A  Monte  Carlo  weighting  technique  developed  by  Schaefer  (1990)  was  utilized  to  determine 
the  most  probable  combination  of  these  parameters.  From  the  results  of  the  Monte  Carlo  simulations, 
high,  medium,  and  low  Bulking  Factors  were  selected  as  a  range  of  probable  values  for  the  sensitivitv 
analysis  that  was  conducted  by  HEC. 

Final  breaching  and  routing  simulations  were  conducted  based  on  what  is  referred  to  throughout  the 
remainder  of  the  report  as  "the  HEC  Breaching  and  Bulking  Scenarios."  They  represent  the  breaching  and 
hulking  characteristics  recommended  and  agreed  upon  by  the  Corps  of  Engineers  and  the  U  S.  Fot.  i 
Service. 


2.1  Physical  Setting 

The  South  Fork  of  Castle  Creek  is  a  perennial  stream  that  drains  an  area  of  2.5  square  miles  on  the 
northwest  flank  of  Mount  St.  Helens.  The  Castle  Lake  blockage  is  located  at  the  confluence  of  South 
Castle  Creek  and  Castle  Creek,  a  tributary  to  the  North  Fork  Toutlc  River,  approximately  48  and  60  miles 
upstream  from  the  communities  of  Castle  Rock  and  Longview- Kelso, W A,  respectively.  Castle  Creek  was 
blocked  by  a  debris  avalanche  that  occurred  during  the  May  18,  19X0  eruption  of  Mount  St.  Helens.  WA. 
Figures  2  and  3  show  the  approximate  pre-eruption  and  post-eruption  topography  of  the  Castle  Lake  studv 
area.  Avalanche  materials  formed  a  blockage  approximately  2,000  lcel  long  at  the  crest  and  is  bounded  bv 
bedrock  ridges  on  either  end  and  averages  about  1,4(X)  feet  wide  from  lake  shore  to  the  downstream  toe. 
Figure  4  shows  a  typical  cross  section  taken  through  the  debris  blockage.  The  blockage  has  a  maximum 
height  of  190  feet  measured  from  the  crest  to  the  toe  and  80  feet  from  the  crest  to  the  lake  surface. 

The  blockage  consists  of  two  major  lithologic  units  which  the  USGS  refers  to  as  the  ancestral  dacite 
unit  and  the  modern  dacitc,  andesite,  and  basalt  unit.  The  modern  (1980  eruption)  materials  are  unsorted, 
mostly  unstratifled  mixtures  of  avalanche  materials,  ranging  in  size  from  sill-  and  clay-sized  particles  to 
large  clasts  more  than  5  feet  in  diameter.  Slopes  from  the  crest  toward  the  lake  are  uniform  and  average 
IV  on  4H.  Slopes  from  the  crest  toward  Castle  Creek  are  more  varied  and  range  from  IV  on  10H  to  IV 
on  2H,  with  the  steepest  downstream  slopes  on  the  western  edge  of  the  blockage.  Vertical  thickness  of  the 
debris  blockage  ranges  from  0  to  250  feet  and  averages  more  than  50  feet  thick.  Location  of  the  deepest 
zone  of  avalanche  materials  corresponds  to  the  former  location  of  the  pre-eruption  South  Fork  Castle 
Creek  alignment.  It  is  believed  that  the  old  South  Fork  Castle  Creek  alignment  resulted  from  breaching 
and  erosion  of  prehistoric  avalanche  deposits  that  formed  there  during  an  earlier  eruption  and  blockage 
sequence  that  occurred  some  2, (XX)  to  3,(XX)  years  ago.  The  original  surface  of  ihe  prehistoric  valley 
deposits  was  eroded  except  for  the  flat  swampy  area  referred  to  as  Castle  Creek  Marsh  in  Figure  2.  The 
marsh  may  have  been  a  remnant  from  the  former  prehistoric  Castle  Lake  bed. 


Figure  2 

Pre-Eruption  Topography  of  the  Study  Area 


Figure  3 

Post  Eruption  Topography  of  the  Study  Area 
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Figure  4 

Generalized  Geologic  Section  Through  the  Castle  Lake  Debris  Blockage 


Following  !he  1980  eruption,  a  new  lake  began  forming  directly  behind  the  debris  blockage  materials 
and  attained  a  volume  of  approximately  19, (XX)  acre-feet  before  an  emergency  spillway  could  be  constructed 
by  the  Corps  of  Engineers  in  1981  to  prevent  possible  overtopping.  Installation  of  the  spillwav  at  the 
eastern  edge  of  the  blockage  stabilized  the  lake  elevation  at  2,577  feet  above  NGVD  (sec  Figure  4).  At 
this  elevation,  the  maximum  depth  in  the  lake  is  110  feet  deep  and  contains  approximately  18,500  acre-feet 
of  water. 


2.2  Characteristics  of  Landslide  Dams 

According  to  Costa  and  Schuster  (1986)  landslide  dams  form  in  a  wide  range  of  physiographic  settings. 
The  most  common  types  of  mass  movements  that  can  form  landslide  dams  include  soil  slumps  and  slides; 
mud,  debris  and  earth  flows;  and  rock  and  debris  avalanches  such  as  those  that  occurred  during  the  1980 
Mount  St.  Helens  eruption.  The  most  common  initiation  mechanisms  for  potential  dam-forming 
landslides  arc  excessive  rainfall,  rapid  snow  melt,  earthquakes  and  volcanic  eruptions. 

Figure  5  shows  that  most  landslide  and  debris  blockage  dams  are  very  short  lived.  Costa  and  Schuster 
(1986)  report  that  for  the  63  documented  cases  they  studied,  22  percent  of  the  landslide  dams  failed  in  less 
than  1  day  after  formation  and  that  half  failed  within  a  period  of  10  days.  Less  than  10  percent  of  the 
natural  debris  blockage  dams  last  more  than  1  year.  They  also  report  that  the  most  frequent  mode  of 
failure  with  debris  blockage  dams  is  by  overtopping.  Figure  6  is  adapted  from  Costa  and  Schuster,  1986 
and  shows  that  more  than  50  percent  of  the  documented  debris  and  landslide  dams  failed  due  to 
overtopping.  The  occurrence  of  a  particular  dam  failure  and  the  magnitude  of  resulting  Hoods  are 
predicated  by  the  siz.e  of  the  blockage,  its  geometric  characteristics  (size  and  depth  of  the  impoundment, 
and  size  and  shape  of  the  blockage),  the  properties  of  the  blockage  materials,  the  rate  of  filling  of  the 
impoundment,  the  volume  of  the  trapped  water,  bedrock  or  engineered  controls  such  as  spillways,  tunnels 
and  diversions. 

The  Castle  Lake  blockage  was  ten  years  old  in  May,  1990  and  appears  to  be  stable  under  its  present 
conditions.  Groundwater  levels  in  the  blockage  and  seeps  along  the  downstream  face  of  the  blockage  have 
been  monitored  since  the  eruption.  According  to  the  Corps’  Geotechnical  Branch  (personal 
communication,  1990),  they  have  seen  no  field  evidence  of  unstable  conditions  in  the  blockage  materials 
since  the  installation  of  the  spillway.  The  Corps  of  Engineers  "Engineering  Analysis  and  Alternative 
Evaluation"  report  (1988)  concludes  that:  (1)  the  risk  associated  with  a  single  event  leading  to  the  failure 
of  Castle  Lake  in  its  existing  condition  is  low;  (2)  the  existing  blockage  is  significantly  larger  than  the 
"minimum  embankment  section"  necessary  to  safely  retain  Castle  Lake;  (3)  local  areas  of  instability  exist 
within  the  blockage,  however,  these  areas  are  outside  the  minimum  embankment  section;  (4)  no  realistic 
failure  mechanism  or  scenario  could  be  developed  that  would  lead  to  the  sudden,  catastrophic  release  of 
Castle  Lake  based  upon  assumed  parameters  (Due  to  the  large  number  of  variables  involved  and 
uncertainties  associated  with  each,  however,  it  is  not  possible  to  completely  eliminate  all  risk),  and  (5)  the 
blockage  exists  in  an  environment  where  rapid  changes  arc  possible  (erosion,  earthquakes,  floods,  volcanic 
eruptions,  etc.).  Monitoring  and  maintenance  arc  necessary  to  ensure  that  the  design  assumptions  that  led 
to  the  above  conclusions  remain  valid. 

The  investigation  reported  herein  is  intended  to  estimate  the  potential  for  flooding  downstream  from 
the  Corps’  Sediment  Retention  Structure  for  various  hypothetical  lake  breaching  scenarios.  It  is  not  the 
intent  of  his  investigation  to  evaluate  any  aspect  of  the  risk  of  failure.  It  merely  quantifies  the  downstream 
flood  potential  for  various  hypothetical  breaching  scenarios. 
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Percentage  of  Dams  Which  Have 
Failed  Below  Indicated  Age 


2J  Information  and  Data  Sources 


Data  collection  for  the  Phase  1  portion  of  this  investigation  began  with  a  field  investigation  of  the 
Castle  Lake  blockage  and  the  downstream  channels  in  November  1989.  The  field  investigation  provided  a 
realistic  perspective  on  the  characteristics  of  the  debris  blockage,  the  location  of  the  historical  South  Fork 
Castle  Creek,  the  characteristics  and  physical  features  of  the  downvalley  deposits  and  the  amount  of 
sediment  and  debris  available  for  llow  bulking  during  high  flow  events.  Many  photographs  were  taken  of 
the  channel  and  valley  sections  along  Castle  Creek  below  the  blockage.  Manning's  n-values  were  estimated 
for  different  sections  along  the  channel  and  overbanks.  Detailed  maps  and  aerial  photographs  were 
obtained  from  the  Portland  District  USACE  and  from  the  U.S.  Forest  Service.  The  maps  and  aerial 
photos  cover  the  area  frem  the  Castle  Lake  debris  blockage  down  to  the  SRS.  Surveyed  cross  sections 
were  also  obtained  from  the  Portland  District  USACE  and  from  the  USGS  Open-File  87-549  by  Meyer 
and  Dodge  (1988).  These  data  were  used  to  develop  47  cross  sections  between  Castle  Lake  and  the  SRS 
and  90  cross  sections  from  the  SRS  to  the  Columbia  River.  The  cross  sections  describe  the  channel  and 
valley  morphology  required  by  the  unsteady  flow  routing  model  developed  for  the  study  reach  from  Castle 
Lake  to  the  SRS  and  from  the  SRS  to  the  Columbia  River.  Dimensions  and  detailed  hydraulic 
information  about  the  Sediment  Retention  Structure  (SRS)  and  its  spillway  and  outlet  works  were  taken 
from  the  Corps’  final  design  manuals  for  the  structure.  These  data  included  volume-elevation  information, 
recently  surveyed  sediment  levels  behind  the  structure,  and  elevation-outflow  information  for  the  low  How 
conduits  and  emergency  spillway. 

During  the  course  of  this  study,  it  was  determined  that  updated  soils  information  was  required  to 
properly  estimate  the  flow  bulking  potential  that  exists  in  debris  deposits  downstream  from  Castle  Lake. 
The  Portland  District  sent  a  soils  investigation  team  to  collect  soil  samples  and  to  measure  in  situ  material 
properties  along  the  channel  and  overbank  areas  below  Castle  Lake.  Data  were  also  collected  along  the 
debris  blockage  itself.  This  information  was  used  to  determine  the  breaching  characteristics  (critical  time 
to  breach  and  ultimate  breach  dimensions)  for  the  debris  blockage.  These  data  were  also  used  to  estimate 
the  range  of  possible  bulking  factors  that  could  occur  downstream  from  the  blockage  due  to  various 
breaching  scenarios. 

Information  and  data  used  during  this  investigation  and  presented  in  this  report  were  obtained  from 
reports,  papers  and  materials  provided  to  HEC  by  the  Portland  District  Corps  of  Engineers,  the  U.S. 

Forest  Service,  the  U.S.  Geological  Survey  and  the  Washington  Stale  Department  of  Ecology.  An 
annotated  bibliography  and  list  of  references  are  provide  in  appendix  A  of  this  report. 

2.4  Breaching  Characteristics  of  the  Debris  Blockage 

Figures  5  and  6  show  that  most  debris  blockage  lakes  fail  within  one  year  of  their  formation  and  the 
most  common  failure  mechanism  is  by  overtopping.  By  installing  an  emergency  spillway  in  1981,  the 
Corps  of  Engineers  essentially  eliminated  the  possibility  of  an  overtopping  failure  of  the  Castle  Lake 
blockage.  Under  present  conditions,  failure  of  the  debris  blockage  would  most  likely  occur  due  to  a 
"piping  type  failure,"  or  as  a  result  of  an  earthquake  occurring  in  conjunction  with  a  severe  hydrologic 
event  that  may  lead  to  a  "heave  type  failure"  (Laenen  and  Orzol,  1987).  One  of  the  first  tasks  of  the  phase 
2  portion  of  this  investigation  was  to  estimate  the  range  of  possible  breach  sizes  and  critical  breach  times. 
Three  different  methods  were  used  to  determine  possible  breach  sizes  and  times.  The  first  two  methods 
arc  statistically  derived  regression  equations,  formulated  by  MacDonald  and  Langridge-Monopolis  (1984) 
and  by  Froclich  (1987).  Both  sets  of  equations  are  based  on  actual  data  from  dozens  of  historic  dam 
failures.  The  MacDonald  and  Langridge-Monopolis  study  was  based  on  data  from  42  man  made  earth  and 
rockfill  dams  (30  carthfill  and  12  that  were  a  combination  of  earth,  clay  cores,  rock  fill,  and  concrete 
faces).  The  Froelich  study  included  data  from  43  man-made  and  landslide  formed  earth  dams.  Both 
studies  resulted  in  a  set  of  graphs  and  equations  that  can  be  used  to  predict  the  approximate  si/c  of  the 
breach  and  the  time  it  takes  for  the  breach  to  reach  its  full  failure  size. 
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The  third  approach  for  estimating  breaching  characteristics  of  the  debris  blockage  was  a  rhvsicaUy  '  im-.I 
computer  mode!  called  BREACH,  developed  by  Dr.  Danny  Frtad  (1989)  of  the  National  Weather  Service. 
The  breach  model  uses  sediment  transport  and  hydraulic  routing  equations  to  simulate  the  formulation  of 
either  a  piping  or  over-topping  type  of  failure.  The  BREACH  computer  model  requires  information  about 
the  physical  dimensions  of  the  dam,  as  well  as  a  very  detailed  description  of  the  soil  properties  of  the  dam 
or  blockage  materials.  Required  soils  information  included: 

1.  D50  (mm) 

2.  Porosity 

3.  Unit  Weight  (lb/ft1) 

4.  Internal  Friction  Angle 

5.  Cohesive  Strength  (lb/ft2) 

6.  D90/D30 

These  parameters  can  be  specified  separately  for  the  inner  core  and  outside  bank  materials  of  a  dam.  In 
the  case  of  the  Castle  lake  blockage,  the  inner  core  material  was  assumed  to  be  the  same  as  the  outer 
banks.  For  this  study,  the  parameters  were  calculated  from  soil  samples  taken  by  the  USGS  and  the 
Portland  District  of  the  Corps  of  Engineers.  A  range  of  appropriate  values  was  extracted  from  the  field 
data.  A  sensitivity  analysis  was  performed  to  see  if  the  BREACH  model  would  predict  different  breach 
sizes  for  different  combinations  of  the  parameters.  The  sensitivity  analysis  showed  that  the  size  of  the 
breach  did  not  vary  significantly  over  the  range  of  parameters  extracted  from  the  field  data.  Tabic  1  shows 
the  range  of  values  used  in  the  sensitivity  analysis,  and  the  final  set  of  values  used  for  this  study. 


Table  1.  Major  soil  properties  used  in  BREACH  model. 


PARAMETER 

RANGE  OF  VALUES 

VALUE  USED 

1.  D50  (mm) 

1.0 -9.0 

1.0 

2.  Porosity 

.34  -  .40 

.38 

3.  Unit  Weight  (lb/ft1) 

100-  145 

125 

4.  Internal  Friction  Angle 

34-36 

35 

5.  Cohesive  Strength  (lb/ft2) 

1  -  400 

2(H) 

6.  D90/D30 

10  -  125 

75 

The  breaching  characteristics  for  each  method,  along  with  the  USGS  heave  scenario  developed  by 
Lacnen  and  Orzol  (1987),  are  summarized  in  Table  2.  Also  shown  are  the  resulting  clear  water  peak  flows 
that  would  occur  for  the  respective  breaching  scenarios.  Costa  and  Schuster  (1988)  developed  a  set  of 
curves  showing  the  potential  energy  of  the  lake  water  versus  peak  discharge  from  historical  dam  failures  of 
various  types  of  dams.  Figure  7  presents  the  five  different  curves  developed  by  Costa  and  Schuster  (198S) 
for  (1)  constructed  dams,  including  earth  and  rockfill  dams,  (2)  landslide  dams,  (3)  Moraine  dams,  (4) 
Glacier  dams,  and  (5)  an  upper  envelope  curve  for  all  of  the  observed  dam  failure  data.  Table  2  lists  the 
peak  discharge  estimated  from  the  Costa  and  Schuster  envelope  curve  (566,000  cfs)  using  the  physical 
characteristics  of  Castle  Lake.  The  U.S.  Bureau  of  Reclamation  (1977)  prepared  an  earlier  curve  of 
observed  peak  discharge  versus  the  hydraulic  depth  of  a  dam  prior  to  failure.  Figure  8  presents  the  U.S. 
Bureau  of  Reclamation’s  curve  and  shows  that  for  an  initial  depth  of  water  behind  a  full  Castle  l.akc,  the 
estimated  peak  discharge  would  be  approximately  370,000  cfs.  The  clear  water  peak  flows  predicted  by  all 
of  six  of  these  different  methods  docs  not  account  for  the  inclusion  of  sediment  front  the  breach. 
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Table  2 


Summary  of  Breaching  Characteristics 


BREACHING 

METHOD 

BOTTOM  WIDTH 
(ET) 

t 

SIDE  SLOPES 
(H/V) 

CRITICAL 
BREACH  TIME 
(HOURS) 

PEAK  11.0 W  EKOM 
CASTLE  LAKE 
(CES) 

U.SX.S  HEAVE1 

675 

1.0 

0.25 

1,510.000 

BREACH  MODEL2 * 
(HEC  SCENARIO) 

480 

0.31 

0.25 

1.180.000 

EROELICH4 

EQUATIONS 

305 

0.31 

0.36 

761.300 

POTENTIAL 
ENERGY4  VERSUS 
PEAK  Q 

RELATIONSHIPS  - 
HISTORICAL  DATA 

- 

- 

- 

566,000 

U.S.B.R  CURVE  - 5 
HISTORICAL  DATA 

- 

- 

- 

370.0OO 

MacDONALD6 

LANGRIDGE- 

MONOPOLIS 

25 

0.31 

0.50 

147,600 

1  lacnen  and  Orzol  (1987) 

2  MI'C's  Breaching  Scenario  using  ihe  NWS  BREACH  model  (1989) 

1  Eroelich  (1987) 

4  Costa  and  Schuster  (1988) 

5  U.S.  Bureau  of  Reclamation  (1977) 

6  MacDonald  and  Langridge-Monopolis  (1984) 


Figure  7 

Potential  Energy  Versus  Peak  Discharge  Relationships  for  Various  Types  of 
Dam  Failures  With  Estimated  Peak  Discharges  for  the  Castle  Lake  Blockage 
(adapted  from  Costa  &  Schuster,  1988) 
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Figure  8 

Envelope  Curve  of  Discharges  Experienced  from  Historical  Dam  Failures 
(from  U.S.  Bureau  of  Reclamation,  1977) 


After  examining  the  physical  characteristics  of  the  Castle  Lake  blockage  and  comparing  the  Head  ig 
characteristics  and  peak  discharge  estimates  from  the  various  breach  development  methods  listed  in  Table 
2,  results  from  the  BREACH  model  are  thought  to  be  the  most  reliable  and  his  method  for  estimating  the 
breaching  characteristics  of  Castle  Lake  the  most  dependable.  'Hie  MacDonald  and  Langridge-Mi  nopolis 
results  arc  considered  inappropriate  because  their  data  set  did  not  include  any  landslide  or  debris  blockage 
failures.  Results  from  the  BREACH  model  were  selected  over  the  Froelich  equations  because  the 
BREACH  model  is  a  more  physically  based  method  and  it  accounts  for  the  material  characteristics  of  the 
blockage  more  explicitly,  therefore  providing  a  better  representation  of  the  specific  problem  at  Castle 
Lake.  Results  from  the  HEC  breaching  scenario  (using  the  BREACH  model)  arc  more  conservative  than 
those  from  other  methods. 

Dozens  of  different  types  of  failure  scenarios  were  tested  to  try  to  duplicate  the  characteristics  of  the 
USGS  proposed  "heave  type"  failure.  No  reasonable  set  of  failure  parameters  could  reproduce  a  breach 
with  a  top  width  of  1000  feet  and  an  approximate  depth  of  175  feet  within  the  critical  breach  time  of  15 
minutes.  HEC  even  tried  to  simulate  the  retrogressive  heave  failure  mechanism  by  initially  removing  65 
percent  of  the  thickness  of  the  blockage  materials  from  the  downstream  face  of  the  blockage,  and  then 
starting  the  BREACH  model  in  a  piping  mode  with  a  full  lake.  Even  under  these  simulated  heave  failure 
conditions,  the  resulting  breach  size  was  similar  in  overall  dimensions  to  those  produced  by  the  "HEC 
breaching  scenario."  After  considerable  discussion  with  engineers  and  geologists  from  the  Corps  of 
Engineers  and  the  U.S.  Forest  Service,  it  was  decided  to  adopt  HEC’s  proposed  breaching  scenario  using 
the  BREACH  model  as  the  most  representative  breaching  approach  for  the  remainder  of  the  investigation. 


2.5  Flow  Bulking  and  Mudflow  Routing  Procedures 

The  BREACH  model  was  used  to  generate  the  dambreak  outflow  hydrographs  from  Castle  Lake  for 
various  initial  lake  levels  and  breaching  scenarios.  Routing  of  the  bulked  dambreak  hydrographs 
(mudflows)  65  miles  from  Castle  Lake  down  the  North  Fork  Toutle  and  Cowlitz  Rivers  to  the  Columbia 
River  was  accomplished  with  Fread’s  (1989)  DAMBRK  model.  Two  separate  routing  reaches  were 
established  as  shown  in  Figure  9.  The  first  reach  (routing  reach  1)  extends  from  Castle  Dike,  16  miles 
down  to  the  SRS.  Flow  bulking  and  dcbulking  processes  occur  within  this  reach  (sec  Lacncn  and  Orzol, 
1987).  The  second  routing  reach  (routing  reach  2)  extends  from  the  SRS  all  the  way  to  the  Columbia 
River  below  Kelso-Longvicw.  The  "mudflow  routing  option"  in  the  DAMBRK  computer  program  was 
used  to  simulate  the  non-Newtonian  hypcrconccntraled  flow  properties  of  the  bulked  discharges  in  routing 
reach  1  downstream  from  Castle  Lake.  The  program  requires  the  user  to  specify  the  expected  mudflow 
properties,  such  as  viscosity  and  initial  shear  strength  of  the  fluid.  Based  on  these  expected  fluid 
properties  and  the  hydraulic  characteristics  of  the  flow,  the  "mudflow  routing  option"  adjusts  the  effective 
friction  loss  terms  in  the  momentum  equations  to  simulate  the  effects  of  hypcrconcentratcd  (bulked)  flow. 

Flow  bulking  is  the  process  whereby  extremely  high  energy  flood  flows  incorporate  additional  bed 
material  (sediment  and  debris)  into  the  flow  by  erosion  and  entrainment,  thus  increasing  the  total  volume 
of  the  flood.  As  the  concentration  of  suspended  material  increases  beyond  a  threshold  of  approximately  20 
to  60  percent  by  volume  (Beverage  and  Culbertson,  1964),  the  fluid-sediment  mixture  begins  to 
demonstrate  non-Newtonian  fluid  characteristics  (i.e.,  the  flow  characteristics  become  dependent  upon  the 
concentration  of  suspended  material).  Many  researchers  have  estimated  that  flow  bulking  can  increase  the 
overall  volume  of  a  hypothetical  dambreak  from  Castle  Lake  from  2  to  5  times  its  original  volume  (Costa, 
1984,  Laencn  and  Orzol,  1989,  and  Scott,  1988).  The  actual  amount  of  flow  bulking  that  occurs  during  an 
event  depends  on  many  factors  and  is  difficult  to  estimate  ahead  of  time.  Schaefer  (1990)  derived  a 
simplified  relationship  that  estimates  the  "ultimate  bulking  factor"  (BF),  given  the  representative  "in  situ" 
soil  characteristics  of  the  avalanche  materials  in  the  valley  and  channels  downstream  front  the  location  of 
hypothetical  lake  breakout.  Soil  porosity  or  void  ratio,  along  with  the  percent  saturation  of  the  soil  are 
required,  along  with  an  estimate  of  the  ultimate  sediment  concentration  by  volume  that  may  occur  during 
the  bulking  process.  Scott  (1985a,  1985b),  Costa  (1984),  Schuster  (1986)  and  Lacncn  and  Orzol  (1987) 
suggest  that  the  suspended  sediment  concentrations  for  this  type  of  an  event  may  range  from  45  to  55 
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percent  by  volume.  Tnereforc,  a  value  of  50  percent  lor  a  representative  ultimate  concentration  '..en  to 
provide  a  reasonable  assumption  for  (his  study 

Bulking  and  debulking  are  likely  to  occur  within  reach  1  from  Castle  Lake  down  to  the  SRS  atcoiding 
to  Laenen  and  Orzol  (1987).  They  also  suggest  that  bulking  will  occur  from  Castle  Lake  to  the  N-l 
structure  and  debulking  (loss  of  suspended  sediment  materials  from  the  flow)  from  the  N-l  to  the  SRS 
(see  Figure  9).  Reach  1  upstream  from  the  N-l  structure,  has  the  deepest  avalanche  deposits,  most 
erodible  channel  materials,  the  narrowest  valley  sections,  the  steepest  stream  slopes,  and  the  greatest 
sediment  transport  potential  Debulking  is  likely  to  occur  downstream  from  the  N-l  structure  where  the 
valley  widens  to  several  thousand  feet,  the  effective  channel  slope  decreases  considerably  and  the  sediment 
transport  capacity  decreases  below  that  necessary  to  sustain  the  high  concentrations  of  materials  entrained 
from  bulking. 

Bulking  of  the  dambreak  flows  up  to  concentrations  of  50  percent  by  volume  as  suggested  by  Laenen 
and  Orzol  (1987)  was  simulated  by  adding  a  series  of  lateral  mudflow  hydrographs  to  the  main  dambreak 
hydrograph  as  it  was  routed  downstream  from  Castle  Lake.  In  this  way,  the  effects  of  flow'  bulking  arc 
essentially  "blended  into  the  main  dambreak  hydrograph"  as  it  moves  down  valley.  This  process  wits 
accomplished  by  first  calculating  the  potential  bulking  factor  for  each  breaching  scenario  using  Schaefer's 
(1990)  method.  Clear  water  dambreak  hydrographs  from  Castle  Lake  for  the  various  breaching  scenarios 
were  then  routed  down  to  the  SRS  with  no  accounting  of  the  potential  bulking  or  debulking  processes. 
Hydrographs  at  several  points  along  bulking  reach  1  (see  Figure  9)  were  extracted  from  the  clear  water 
runs.  These  hydrographs  were  then  multiplied  by  a  ratio  (the  estimated  subreach  Bulking  Factor  (BF)]  to 
account  for  the  bulked  volume  of  material  that  would  need  to  be  added  to  each  subreach  to  obtain  me 
overall  flow  bulking  and  total  flow  volume  for  the  entire  bulking  reach  for  each  breaching  scenario. 
Distribution  of  the  amount  of  bulked  material  entering  the  flow  with  distance  in  reach  1  was  allocated 
according  to  the  longitudinal  distribution  of  sediment  transport  capacity  along  the  bulking  reach.  This 
method  produces  a  more  realistic  (nonlinear)  relationship  between  local  hydraulic  conditions  (depth,  width 
and  velocity)  and  the  longitudinal  change  in  the  amount  of  bulked  flow  to  be  blended  into  each  subreach. 

The  tendency  for  flow  bulking  actually  increases  for  a  short  distance  downstream  from  the  blockage 
because  the  valley  is  relatively  narrow  and  very  steep  and  the  transport  capacity  is  very  high  along  the  front 
of  the  dam  break  bore.  Consequently,  as  the  flow  picks  up  more  material  from  bulking  processes,  the 
effective  discharge  also  goes  up  until  the  valley  widens  and  flattens  enough  to  begin  attenuating  the  flow. 
This  phenomenon  of  the  magnitude  of  the  flow  increasing  due  to  bulking  for  some  distance  downstream 
from  the  initial  breakout  is  clearly  demonstrated  in  the  results  shown  and  discussed  in  later  sections  of  this 
report.  As  the  valley  widens  and  the  channel  slope  decreases,  the  amount  of  bulking  goes  down  in 
proportion  to  the  reduced  transport  capacity  for  that  subreach.  This  approach  was  used  to  distribute  the 
amount  of  bulked  material  (nonlincarly)  into  the  flows  from  the  dam  break  hvdrograph  as  it  was  routed 
dynamically  downstream  from  Castle  Lake  to  the  N-l  structure.  The  last  step  was  to  run  the  model  again 
using  the  "mudflow  option"  with  the  estimated  mudflow  properties  for  each  reach,  while  adding  (blending) 
the  lateral  bulked  flow  hydrographs  to  main  flow  to  account  for  the  dynamic  flow  bulking. 

Laenen  and  Orzol,  (1987)  suggest  that  debulking  occurs  in  lower  portion  of  reach  l,  from  the  N-l 
structure  to  the  SRS.  The  valley  widens  rapidly  along  this  part  of  the  reach  to  several  thousand  feet  wide, 
the  stream  slope  decreases,  and  there  is  insufficient  sediment  transport  capacity  to  sustain  continued 
bulking  of  the  flows.  The  debulking  process  is  handled  in  a  similar  manner  as  was  the  bulking  process. 
Hydrographs  from  the  Clearwater  runs  for  the  different  breaching  scenarios  are  used  to  establish  the  shape 
and  timing  of  lateral  flows  that  would  be  extracted  (dcbulked)  from  the  main  flood  wave.  Once  again  the 
distribution  of  the  amount  of  debulking  is  allocated  along  (he  debulking  reach  in  proportion  to  the  change 
in  the  amount  of  sediment  transport  capacity  occurring  along  the  reach.  Debulking  is  simulated  with  the 
DAMBRK  model  by  using  negative  lateral  flow  hydrographs  to  remove  the  amount  of  flow  volume  lost 
due  to  debulking.  Laenen  and  Orzol  (1987)  used  similar  methods  to  simulate  the  expected  debulking 
below  the  N-l  structure.  This  process  was  repeated  for  each  breaching  scenario  that  was  evaluated. 
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2.6  Estimating  the  Ultimate  Bulking  Factor  for  Various  Dambreak  Scenarios 


The  magnitude  of  the  ultimate  bulking  factor  is  not  only  a  function  of  the  characteristics  of  the 
dambreak  hydrograph  producing  the  flow,  but  also  a  function  of  the  in  situ  properties  of  the  valley  and 
channel  deposits  where  the  dambreak  flows  will  occur.  Examination  of  field  data  collected  by  the  Corps  of 
Engineers  and  the  USGS,  shows  that  actual  field  values  for  porosity  and  void  ratio,  percent  saturation  and 
ultimate  suspended  sediment  concentration  can  vary  according  to  the  season,  material  type  and  location 
along  the  channel.  For  the  purposes  of  this  investigation,  a  method  was  derived  to  bracket  the  range  of 
possible  material  properties  observed  in  the  field  and  to  assign  a  level  of  confidence  to  the  many  possible 
combinations  of  the  three  main  variables  (porosity,  percent  saturation  and  ultimate  suspended  sediment 
concentration)  that  can  occur.  For  this  study,  the  following  range  of  values  for  the  main  parameters  was 
agreed  upon: 

1.  Porosity  (0.25  -  0.45) 

2.  %  Saturation  (45  -  90) 

3.  Ultimate  concentration  (0.30  -  0.60) 

Schaefer  (1990)  developed  a  statistical  weighting  method  using  Monte  Carlo  sampling  techniques,  to 
estimate  the  magnitude  of  the  "ultimate  bulking  factor"  according  to  observed  ranges  in  the  magnitudes  of 
the  three  main  variables  used  to  compute  the  bulking  factor  (BF).  The  Monte  Carlo  method  computed  a 
"cumulative  probability"  of  99.7%  to  the  USGS’  "heave  scenario  with  a  bulking  factor  of  4.46."  The 
method  assigns  a  95%  cumulative  probability  to  the  HEC  "piping  scenario  with  a  BF  =  3.32,"  and  50'  <  o 
the  HEC  "piping  scenario  with  a  BF  =  2.5."  The  cumulative  probability  means  that  of  all  the  possible 
bulked  flows  that  can  occur  during  a  breaching  of  Castle  Lake  with  the  range  of  observed  material 
properties  in  the  channel  downstream  from  the  blockage,  99.7,  95,  and  50  percent  of  the  events  will  have 
bulking  factors  less  than  4.46,  3.32  and  2.5,  respectively.  Therefore,  of  all  the  different  combinations  of 
values  of  the  three  bulking  factor  variables,  less  than  1/2  of  a  percent  of  all  the  possible  combinations 
would  yield  a  bulking  factor  as  large  as  4.46  for  the  USGS’  breaching  scenario.  Approximately  95  percent 
of  all  the  possible  events  would  be  less  than  the  3.32  bulking  factor  associated  with  HEC’s  piping  scenario 
and  about  50  percent  of  all  the  events  would  have  a  bulking  factor  of  2.5  or  less.  This  range  of  bulking 
factors  was  used  to  evaluate  the  potential  for  flooding  downstream  from  the  SRS  for  three  hypothetical 
lake  breaching  scenarios.  Results  and  conclusions  from  the  analyses  are  presented  in  following  sections  of 
the  report. 


2.7  Other  Boundary  Conditions  Considered 

Figure  9  presents  a  schematic  diagram  of  the  approximate  location  and  extent  of  the  routing  reaches 
considered  during  this  investigation.  Three  different  breaching  scenarios  were  considered  for  the 
hypothetical  failure  of  the  Castle  Lake  blockage.  Tables  2  and  3  summarize  the  specific  breaching 
scenarios  that  were  considered  along  with  the  boundary  conditions  and  range  of  bulking  factors  that  were 
applied.  Tables  2  and  3  also  list  the  key  locations  along  routing  reaches  1  and  2  where  computed  results 
are  displayed  for  comparison  and  discussion.  As  discussed  in  section  2.5,  routing  of  the  hypothetical 
breaching  events  was  carried  out  in  two  different  reaches.  Reach  1  included  the  flow  bulking  and 
dcbulking  effects  on  the  flows  from  Castle  Lake  to  the  SRS.  Reach  2  extends  from  the  SRS  all  the  way  to 
the  Columbia  River.  It  is  assumed  that  most  of  the  sediment  and  debris  load  will  be  captured  behind  the 
SRS  and  that  flow  bulking  and/or  dcbulking  is  insignificant  throughout  reach  2.  Tables  3  and  4  show  the 
initial  Castle  Lake  water  surface  elevations  that  were  considered:  (1)  lake  full  at  elevation  2,580  NGVD, 

(2)  lake  lowered  30  feet  to  elevation  2,550  NGVD,  and  (3)  lake  lowered  60  feet  to  elevation  2.520  NGVD. 
The  spillway  crest  elevation  at  Castle  Lake  is  2,577  feet  NGVD.  The  starting  lake  elevation  of  2,580 
corresponds  to  the  estimated  lake  level  that  would  occur  during  a  500  year  rainfall  event  positioned  over 
the  3.0  square  mile  drainage  area  above  the  lake.  Under  these  conditions,  a  discharge  of  approximately 
6, (KM)  cfs  would  be  flowing  out  of  Castle  Creek  as  base  flow  prior  to  the  d3m  break.  The  estimated  runoff 
and  base  flow  in  the  North  Fork  Toutle  River  upstream  from  the  N-l  structure  prior  to  the  dam  break  was 
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rABLE  4 


Summary  of  Simulations  Performed  and 
Locations  Where  Results  Were  Printed 


CASTLE  LAKE  TO  SRS  ROUTING  REACH 


USGS  HEAVE  SCENARIO 
LAKE  LEVEL 


BF=4.46 

KEY  LOCATIONS  WHERE 
RESULTS  WERE  PRINTED 


FULL 

-30’ 

X 

0.0.  1.87.  5.31.  10.19.  12.8.  15.71 

X 

0.0,  1.87.  5.31,  10.19.  12.8,  15.71 

SRS  TO 

COLUMBIA 

RIVER 

ROUTING  REACH 

USGS 

HEAVE  SCENARIO 

BF=4.46 

LAKE 

LEVEL 

SRS 

KEY  LOCATIONS  WHERE 
RESULTS  WERE  PRINTED 

FULL 

- 

-30’ 

FULL 

EXIST. 

COND. 

(MILES  D/S  FROM  CAST.  L.) 

X 

X 

15.91,  23.21,  33.91,  47.91.  59.81 

X 

X 

15.91,  23.21,  33.91.  47.91.  59.81 

X 

X 

15.91,  23.21,  33.91,  47.91.  99.81 

X 

X 

15.91.  23  21.  33.91.  47.91.  59.81 
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approximately  20, (XXI  cfs.  Other  assumed  base  flow  conditions  include  (1)  0,(XX)  els  from  the  SRS,  <1 
5, 70()  cfs  from  the  Green  River  into  the  North  Fork  Toutle  River,  (3)  6,100  cfs  from  the  South  Fork 
Toulle  River,  and  (4)  52,0(X)  efs  in  the  Cowlitz  river  just  upstream  from  the  confluence  with  the  Toutle 
River.  These  base  flow  conditions  were  recommended  by  the  Portland  District’s  Hydrology  and  River 
Engineering  Branch  for  the  purposes  of  this  investigation.  The  total  base  flow  in  the  Cowlitz  river 
downstream  from  the  confluence  with  (he  Toutle  River  as  a  result  of  the  cumulative  base  flows  from  the 
rivers  entering  above  is  approximately  69,8(X)  cfs. 

The  other  key  boundary'  condition  considered  as  a  variable  during  the  investigation  was  the  initial 
storage  condition  behind  the  SRS.  If  an  event  were  to  occur  today  (in  September,  1990)  the  storage 
conditions  would  be  represented  by  the  "existing  conditions"  in  the  SRS  (see  Tables  3  and  4  for  routing 
reach  2.  SRS  to  the  Columbia).  However,  the  SRS  is  designed  to  trap  sediment  materials  from  Mount  St. 
Helens  and  the  drainages  upstream  from  the  SRS.  Therefore,  over  time  the  amount  of  available  storage 
will  decrease  until  the  SRS  is  completely  fuil  of  sediment  and  debris  up  to  the  crest  of  the  spillway.  Even 
under  "completely  full  conditions",  the  SRS  provides  a  great  deal  of  storage  and  attenuation  of  flood  waves 
as  will  be  discussed  in  the  next  section. 


3.  Results  and  Discussion 

Figure  10  presents  curves  of  the  computed  peak  discharge  versus  river  mile  resulting  from  the 
hypothetical  failure  of  the  Castle  Lake  blockage  due  to  piping  for  three  initial  lake  levels  and  two  diTercnt 
SRS  conditions.  The  "ultimate  bulking  factor"  for  the  scenarios  simulated  and  presented  in  Figure  It)  was 
2.50.  Results  for  this  scenario  are  considered  to  be  the  approximate  average  (representative  of  50  'r  of 
the  events  that  may  occur)  based  on  the  field  data  presently  available.  Figure  11  presents  similar  curves  of 
peak  discharge  versus  river  mile  for  the  same  type  of  piping  failure  but  with  an  "ultimate  bulking  factor"  of 
3.32.  These  results  represent  a  scenario  with  a  mode  of  failure  and  a  magnitude  of  bulking  that  could  be 
equalled  or  exceeded  by  only  5  percent  of  the  possible  flow  events  that  were  considered  (or  95  %  of  the 
events  will  be  less  than  this  magnitude).  This  breaching  and  bulking  scenario  is  recommended  by  the 
Corps  of  Engineers,  the  U.S.  Forest  Service  and  the  Washington  Stale  Department  of  Ecology  for  the 
evaluating  the  flooding  effects  a  hypothetical  beaching  of  the  Castle  Lake  blockage.  Figure  1 2  presents 
similar  results  for  the  USGS’  "heave  type  breaching  scenario"  introduced  by  Laencn  and  Orzol  (1987)  and 
modified  by  HEC.  For  this  scenario  HEC  used  Laenen  and  Orzol’s  breach  hydrograph,  but  applied  HEC\ 
energy  related  bulking  approach  along  with  Schaefer’s  "ultimate  bulking  factor"  of  4.46  based  on  the 
assumed  field  conditions.  Only  two  initial  lake  elevations  were  considered  for  this  scenario.  The  results 
arc  presented  in  Figure  12  for  lake  full  conditions  and  the  lake  lowered  30  feet.  The  "ultimate  bulking 
factor"  associated  with  the  USGS’  breaching  scenario  was  4.46.  This  represents  a  cumulative  probability 
(sec  Section  2.6)  of  99.7  percent:  c.g.  99.7  percent  of  the  events  considered  will  be  less  than  the  computed 
discharges  presented  in  Figure  12.  Table  5  summarizes  the  routing  results  at  key  locations  below  Castle 
Lake  for  the  three  different  breaching  and  bulking  scenarios  considered:  (1)  HEC’s  piping  failure  with  the 
median  BF  =  2.5,  (2)  HEC’s  piping  failure  with  the  95  %  BF  =  3.32  and  (3)  the  USGS’  heave  failure  w  ith 
the  99.7  %  BF  =  4.46.  The  peak  discharges  listed  in  Table  5  should  approximately  equal  the  peak 
discharges  plotted  in  Figures  10  through  12.  Note  that  the  initial  elevation  of  Castle  Lake  and  amount  of 
sediment  stored  in  the  SRS  make  a  considerable  difference  in  the  magnitude  of  the  initial  dambreak  flow 
and,  therefore,  a  great  deal  of  difference  by  the  time  the  flood  wave  is  routed  to  the  SRS  and  downstream 
to  Castle  Rock  and  Kelso  -  Longview.  For  example.  Table  5  lists  the  routing  results  for  the  Corps’- 
recommended  scenario  (runs  7  through  12:  HEC’s  piping  scenario  with  the  95  %  BF  =  3.32).  Note  that 
with  Castle  Lake  full  (elevation  =  2580  NGVD)  and  existing  at  the  SRS  (see  run  7  in  Table  5),  the  peak 
discharge  in  the  North  Fork  Toutle  River  below  the  SRS  is  reduced  85  percent  (from  695.0(X)  to  105, 200 
cfs)  wi.h  the  SRS  in  place.  If  the  lake  is  lowered  30  feet  or  60  feet  (runs  8  and  9),  the  flow  reduction  due 
to  storage  in  the  SRS  is  approximately  82  and  95  percent,  respectively.  Even  if  the  SRS  were  initially  full 
of  sediment  (see  runs  10,  11  and  12  in  Table  5),  the  peak  discharge  from  a  failure  of  Castle  Lake  would  be 
significantly  reduced.  For  this  case  the  flows  would  be  reduced  by  76,  66  and  56  perccn',  for  the  lake  full, 
the  lake  lowered  30  feet  and  the  lake  lowered  60  feet,  respectively. 
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Figure  10 

Peak  Discharge  Versus  River  Mile  Resulting  from  the 
Hypothetical  Failure  of  Castle  Lake  for  Three  Initial 
Lake  Levels  and  Two  Different  SRS  Conditions 
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Figure  1 1 

Peak  Discharge  Versus  River  Mile  Resulting  from  the 
Hypothetical  Failure  of  Castle  Lake  for  Three  Initial 
Lake  Levels  and  Two  Different  SRS  Conditions 
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Figure  12 

Peak  Discharge  Versus  River  Mile  Resulting  from  the 
Hypothetical  Failure  of  Castle  Lake  for  Two  Initial 
Lake  Levels  and  Two  Different  SRS  Conditions 
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Figures  10  through  12  clearly  show  the  boundary  condition  effects  on  the  peak  discharges  for  vnnr.is 
combinations  of  initial  lake  levels,  SRS  conditions  and  bulking  factor.  Because  the  SRS  substantially 
reduces  the  flow  downstream  from  the  SRS,  dambreak  discharges  added  to  the  base  flows  in  the  Toutlc 
and  Cowlitz  Rivers  downstream  front  the  SRS  may  not  be  significantly  greater  than  the  initially  assumed 
base  flow  conditions.  For  the  case  of  runs  7  through  9  (see  Figure  1 1  and  Table  5).  only  39,600,  28,700 
and  400  cfs  are  added  to  the  base  flow  in  the  Cowlitz  River  at  Kelso  -  Longview  due  to  the  dambreak 
discharges  estimated  with  the  initial  lake  full,  lowered  30  feet  and  lowered  60  feci,  respectively.  If  the  SRS 
were  Ml  of  sediment  at  the  time  of  failuie  (runs  10  liuough  12),  59,000,  46,000  and  18,600  cfs  would  Ke 
added  to  the  69,800  cfs  base  flow  in  the  Cowlitz  River  at  Kelso  -  Longview.  This  represents  an  increase  in 
discharge  above  the  assumed  base  flow  in  the  Cowlitz  River  of  85,  66  and  27  percent,  respectively  for  the 
three  different  initial  lake  conditions. 

Figure  13  summarizes  the  estimated  Hood  magnitudes  from  Castle  Lake  to  Kelso  -  Longview  for  the 
HEC  breaching  and  bulking  scenario  (BF  =  3.32),  for  two  initial  lake  conditions  (full  and  lowered  30  feet) 
and  existing  (1990)  conditions  in  the  SRS.  The  hypothetical  HEC  mudflow  conditions  arc  compared  to 
the  1980  Mount  St.  Helens  mudflow  and  the  probable  maximum  flood  (assuming  the  SRS  was  not  there 
during  the  PMF).  The  spillway  capacity  at  the  SRS  is  indicated  in  Figure  13,  along  with  the  approximate 
maximum  channel  capacities  near  the  communities  of  Castle  Rock  and  Kelso  -  Longview.  Both 
hypothetical  HEC  dambreak  floods  (for  lake  full  and  lowered  30  feet)  represent  "existing  conditions"  at  the 
SRS  and  in  the  Cowlitz  River.  The  routed  dambreak  flows  for  these  "existing  conditions"  are  all  contained 
within  the  present  Cowlitz  River  at  Castle  Rock  and  Kelso  -  Longview.  If  the  SRS  were  full  of  sediment 
initially,  the  estimated  HEC  dambreak  flood  would  exceed  the  present  channel  capacity  near  Castle  Rock 
by  28,700  and  11,800  cfs  if  Castle  Lake  were  initially  full  or  lowered  by  30  feet,  respectively.  The  channel 
near  Kelso  -  Longview  (estimated  to  be  140,000  cfs)  presently  has  sufficient  capacity  to  contain  routed 
dambreak  flows,  even  if  the  SRS  were  full  of  sediment.  Estimated  peak  flows  from  the  lessor  HEC 
scenario  (with  a  bulking  factor  of  2.5)  would  all  be  contained  within  the  present  channel  at  Castle  Rock 
and  Kelso  -  Longview.  If  an  event  as  rare  as  the  hypothetical  "heave  failure  scenario"  (with  a  bulking 
factor  of  4.46)  were  to  happen  under  the  "existing  conditions",  the  channel  capacity  near  Castle  Rock 
would  be  exceeded  for  both  initial  lake  levels  that  were  considered  (Castle  Lake  full  or  lowered  by  30 
feet).  Flood  flows  would  still  be  contained  in  the  present  channel  near  Kelso  -  Longview  under  "existing 
conditions."  If  the  SRS  were  initially  full,  both  Castle  Rock  and  Kelso  -  Longview  could  experience 
flooding  for  lake  full  initial  conditions.  If  the  lake  were  lowered  30  feet,  Castle  Rock  could  still  experience 
flooding,  while  the  flow  at  Kelso  -  Longview  would  be  barely  contained  in  the  channel. 

Therefore,  with  the  present  "existing  conditions"  at  the  SRS  and  in  the  Cowlitz  River,  all  of  the  HEC  - 
recommended  flooding  scenarios  would  be  fully  contained  at  Castle  Rock  and  Kelso  -  Longview.  The 
resulting  flows  would  be  similar  to  a  100  year  flood  event  in  the  Cowlitz  River.  If  the  SRS  were  full  of 
sediment,  all  of  the  HEC  -  recommended  flooding  scenarios  would  be  fully  contained  at  Kelso  -  Longview, 
but  not  at  Castle  Rock  for  cither  lake  full  or  lake  lowered  30  feet  conditions.  None  of  the  hypothesized 
breaching  and  bulking  scenarios  will  exceed  or  overtop  the  SRS  for  either  "existing  conditions"  or  "full 
conditions."  It  is  expected,  however,  that  for  any  major  flooding  event  with  the  magnitude  of  the 
hypothesized  breaching  of  Castle  Lake,  that  significant  quantities  of  sediment  and  debris  will  enter  the 
SRS,  thus  reducing  its  present  storage  capacity  and  active  life. 

During  large  flood  events  the  primary  concern  is  usually  for  channel  capacity  and  whether  the  peak 
discharge  will  be  contained  within  the  existing  channel.  The  routing  results  presented  in  figures  10 
through  13  and  Table  5  show  the  beneficial  effects  of  the  SRS  in  reducing  the  peak  discharges  in  the 
channels  downstream  from  the  SRS.  Table  6  presents  the  average  channel  velocities  computed  at  key 
locations  below  Castle  Lake  for  the  breaching  and  bulking  scenarios  that  were  considered.  Maximum 
velocities  occur  just  below  the  breach.  For  the  HEC  recommended  scenario  (with  a  bulking  factor  of  3.32) 
the  maximum  velocity  just  downstream  from  the  breach  is  approximately  27  feet  per  second.  Velocities  ol 
this  magnitude  occurring  over  the  loose  debris  avalanche  deposits  can  readily  lead  to  the  kinds  of  flow 
bulking  scenarios  described  in  Section  2.5.  By  me  time  the  dambreak  bore  reaches  the  N-l  structure,  it 
has  decreased  its  velocity  by  almost  half  to  14.6  fps.  This  supports  the  debulking  concept  also  presented  in 
Section  2.5.  Downstream  trom  the  N-l,  the  flow  continues  to  slow  down,  but  only  slightly  until  it  enters 
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Table  5 


Routing  Results  at  Key  Locations  Below  Castle 
Lake  for  the  HEC  and  USGS  Breaching  Scenarios 


Initial 

Qc  Just 

Outflow 

Qp  at 

Qp 31 

Initial 

Castle 

Below 

Qp  at  N-l 

Inflow  to 

from 

Castle 

Ixmgvicw 

Run 

SRS 

l^ake 

Castle 

Qp  at 

Structure 

SRS 

SRS 

Rock 

Kelso 

No. 

Condition 

F.levation 

Lake 

RM  5.31 

RM  10.21 

RM  161 

Spillway 

RM  48.01 

RM  59.8 1 

(NGVD) 

(cfs) 

(cfs) 

(cfs) 

(cfs) 

(cfs) 

(cfs) 

(cfs) 

gg  §§§  I 


-X>\J  *'*£<'-  \ 


Peak  Flows,  Qp  (cfs)  for  I  lEC’s 

i  Dam  Breach  (Piping)  Scenario,  with  BF 

=  2.50: 

1 

lixisling 

2,580  (full) 

1,321,300 

1,360,400 

876,900 

351.400 

47,600 

99,800 

82.(XX) 

2 

Existing 

2,550  (-301) 

900,500 

746,400 

523,400 

186,500 

18,800 

77.100 

72.100 

3 

Iucisting 

2,520  (-60’) 

528,900 

350,900 

232,700 

87,100 

6,000 

70.200 

70.2(h) 

4 

Full 

2,580  (full) 

1,321,300 

1,360,400 

876,900 

351,400 

108,900 

143,000 

1 04,<>00 

5 

Full 

2,550  (-30’) 

900,500 

746,400 

523,400 

186,500 

79,500 

122,300 

93,500 

6 

Full 

2,520  (-60’) 

528,900 

350,900 

232,700 

87,100 

40.600 

96,900 

81  ,‘8)0 

Peak  Flows,  Qp  (cfs)  for  HEC’s 

Dam  Breach  (Piping)  Scenario,  with  BF= 

=  3  32’1 

7 

Existing 

2,580  (full) 

1,321,300 

1,540,000 

1,136,000 

695,000 

105,200 

139,700 

109,400 

8 

Existing 

2,550  (-30  ) 

900,500 

990,100 

713,100 

352.500 

62,000 

113.700 

98,500 

9 

Existing 

2,520  (-60  ) 

528,900 

461,400 

340.500 

131,800 

6,000 

70,200 

70.2(H) 

10 

Full 

2,580  (full) 

1,321,300 

1,540,000 

1,136,000 

695,000 

167,200 

178,700 

128,800 

II 

Full 

2,550  (-30’) 

900,500 

990,100 

713,100 

352,500 

119,900 

151,800 

115.800 

12 

Full 

2,520  (-60  ) 

528,900 

461,400 

340,500 

131,800 

57,600 

107,900 

88.400 

Peak  Flows,  Qp  (cfs)  for  i 

the  USGS’s  Heave  Scenario,  with  BF  =4.4b4 

13 

Existing 

2^80  (full) 

1,692,000 

2,311,400 

1,827,900 

1,024,300 

184.8IXJ 

189,900 

134.000 

14 

Existing 

2,550  (-30’) 

1,181,800 

1,642,400 

1,113,700 

708.600 

119,600 

155,000 

113.7(H) 

15 

Full 

2,580  (full) 

1,692,000 

2,311,400 

1,827,900 

1,024.300 

237.000 

2.10,300 

154.8(H) 

16 

Full 

2,550  (-30  ) 

1,181,800 

1,642,400 

1,113,700 

708,600 

179,100 

192..300 

135.1(H) 

’  River  Mile  (RM)  locations  designate  distance  downstream  front  Castle  Lake  in  miles. 

1  llf-C’s  Hulking  Scenario:  BF  =  2.50;  n  =  0.25,  %Sat  =  65%,  Max  Bulk  Cone  =  50%  by  Vol 
1  I  [EC’s  Bulking  Scenario:  BF  =  3.32;  n  =  0.38,  %Sat  =  65%,  Max  Bulk  Cone  =  50%  by  Vol 
1  tISGS’s  Bulking  Scenario:  BF  =  4.46;  n  =  0.38,  %Sat  =  90%,  Max  Bulk  Cone  =  50%  by  Vol 
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Table  6 


Computed  Channel  Velocities  at  Key  Locations  Below 
Castle  Lake  for  the  HEC  and  USGS  Breaching  Scenarios 
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(NGVD) 

Vp  Just 
Below 
Castle 
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I  A>ngview 
Kelso 

RM  59.8 1 
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#1111111111 

\Y‘ 
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Illllill 
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Max  Vels,  Vp  (fps)  for  HEC’: 

s  Dam  Breach  (Piping)  Scenario,  with  BP' 

1!  : 

f  J 

9. 

l 

Existing 

2,580  (full) 

26.13 

20.36 

12.99 

9.30 

7.24 

4.67 

1.44 

2 

Existing 

2,550  (-30  ) 

22.56 

17.13 

10.52 

6.97 

4.98 

3.98 

3.70 

3 

Existing 

2,520  (-60’) 

18.05 

13.87 

7.72 

4.98 

3.13 

3.74 

3.64 

4 

Full 

2,580  (full) 

26.13 

20.36 

12.99 

930 

9.72 

5  68 

"  30 

5 

Full 

2,550  (-30’) 

22.56 

17.13 

10.52 

6.97 

8.74 

5.20 

4.89 

6 

Full 

2,520  (-60’) 

18.05 

13.87 

7.72 

4.98 

6.81 

4.55 

4.43 

Max  Vels,  Vp  (fps)  for  HEC’s 

Dam  Breach  (Piping)  Scenario,  with  Bl'= 

•3.32'* 

7 

Existing 

2,580  (full) 

26.64 

21.65 

14.62 

12.28 

9.61 

5.60 

5.39 

8 

Existing 

2,550  (-30’) 

22.61 

18.72 

12.13 

9.31 

7.99 

4.95 

4.43 

9 

Existing 

2,520  (-60  ) 

18.05 

15.01 

8.83 

5.94 

3.13 

3.74 

3.64 

10 

Full 

2,580  (full) 

26.64 

21.65 

14.62 

12.28 

11.14 

6.52 

609 

11 

Pull 

2,550  (-30’) 

22.61 

18.72 

12.13 

9.31 

10.03 

5.88 

5.64 

12 

Full 

2,520  (-60’) 

18.05 

15.01 

8.83 

5.94 

7.78 

4.83 

4.68 

Max  Vels,  Vp  (fps)  for  the  USGS's  Heave  Scenario,  with  BF=4.464 


13 

Existing 

2,580  (full) 

28.63 

23.61 

16.80 

13.9.3 

11.48 

6.73 

6.27 

14 

Existing 

2,550  (-30’) 

24.91 

22.12 

14.62 

12.37 

10.01 

596 

5.02 

15 

Full 

2.580  (full) 

28.63 

23.61 

16.80 

13.93 

12.38 

7.42 

6.94 

16 

Full 

2,550  (-30’) 

24.91 

22.12 

14.62 

12.37 

11.37 

6.73 

6.33 

1  River  Mile  (RM)  locations  designate  distance  downstream  from  Castle  Lake  in  miles. 

2  lire's  Bulking  Scenario:  BP  =  2.50;  n  =  0.25,  %Sal  =  65%,  Max  Bulk  Cone  =  50%  by  Vol 
!  lire's  Bulking  Scenario:  BF  =  .132;  n  =  0.38,  %Sat  =  65%,  Max  Bulk  Cone  =  50%  by  Vol 

’  USGS’s  Bulking  Scenario:  BP'  =  4.46;  n  =  0.38.  %Sat  =  90%,  Max  Bulk  Cone  =  50%  by  Vol 
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Estimated  Flood  Magnitudes  from  Castle 
Lake  to  Kelso-Longview  for  Present 
Basin  Conditions  (1990) 


the  SRS  storage  pool  area  (velocity  here  is  approximately  12  fps).  Flow  leaving  the  SRS  is  controlled  h\ 
the  spillway  and  outlet  works.  The  maximum  average  channel  velocity  just  downstream  from  the  SRS  is 
approximately  9.6  fps.  By  the  lime  the  Hood  wave  reaches  the  Cowlitz  River  and  mixes  with  the 
cumulative  base  flow  in  the  Cowlitz  (69,800  cfs),  it  slows  to  about  4  to  6  feet  per  second. 

The  SRS  protects  communities  and  the  river  sections  downstream  from  it  but  docs  not  affect  those 
areas  upstream  from  the  SRS.  Therefore,  there  is  no  protection  above  the  SRS  from  a  hypothetical  failure 
of  Castle  Lake.  Because  of  this,  and  the  possibility  that  a  failure  of  Castle  Lake  would  greatly  reduce  the 
sediment  storage  capacity  and  effective  life  of  the  SRS,  additional  field  investigations  pertaining  to  the 
geotechnical  stability  of  the  Castle  Lake  blockage  and  regularly  scheduled  monitoring  are  recommended. 


4.  Conclusions 

The  following  conclusions  arc  drawn  from  the  results  of  this  investigation: 

1)  Numerical  methods  developed  by  the  National  Weather  Service  (Frcad,  1989),  the  Hydrologic 
Engineering  Center,  and  Schaefer  (1990),  were  successfully  used  to  simulate  the  hypothetical 
breaching,  flow  bulking  and  unsteady  mudflow  routing  that  would  result  if  the  Castle  Lake 
blockage  dam  were  to  fail. 

2)  Even  though  debris  blockage  dams  form  in  a  wide  variety  of  physiographic  settings,  most 
debris  blockage  dams  arc  very  short  lived.  Costa  and  Schuster  (1986)  report  that  for  the  63 
documented  cases  they  studied,  22  percent  of  the  landslide  dams  failed  in  less  than  1  day  after 
formation  and  that  half  failed  within  a  period  of  10  days.  Less  than  10  percent  of  the  natural 
debris  blockage  dams  last  more  than  1  year. 

3)  More  than  50  percent  of  the  documented  debris  and  landslide  dams  failed  due  to 
overtopping.  The  occurrence  of  a  particular  dam  failure  and  the  magnitude  of  resulting 
floods  arc  predicated  by:  the  size  of  the  blockage;  its  geometric  characteristics  (size  and  depth 
of  the  impoundment,  and  size  and  shape  of  the  blockage);  the  properties  of  the  blockage 
materials;  the  rate  of  filling  of  the  impoundment;  the  volume  of  the  trapped  water;  and 
bedrock  or  engineered  controls  such  as  spillways,  tunnels  and  diversions. 

4)  The  Castle  Lake  blockage  was  ten  years  old  in  May,  1990  and  appears  to  be  stable  under  its 
past  and  present  conditions.  The  Portland  District  Corps  of  Engineers  installed  an 
emergency  spillway  in  October  of  1981  to  stabilize  the  lake  elevation  at  2577  feet  NGVD. 
Groundwater  levels  in  the  blockage  and  seeps  along  the  downstream  face  of  the  blockage 
have  been  monitored  since  the  eruption.  According  to  the  Corps’  Geotechnical  Branch 
(personal  communication,  1990),  they  have  seen  no  field  evidence  of  unstable  conditions  in 
the  blockage  materials  since  the  installation  of  the  spillway. 

5)  The  estimated  peak  discharge  from  a  hypothetical  failure  of  the  Castle  Lake  blockage  exceeds 
the  peak  discharges  predicted  from  potential  energy  versus  peak  discharge  relationships 
developed  from  historical  dam  failures  by  more  than  2.3  times  (see  Figure  7).  It  exceeds  the 
predicted  peak  discharge  envelope  curve  from  historical  dam  failures  by  3.6  times  (see  Figure 
8).  Therefore,  the  Corps  recommended  breaching  and  bulking  scenario  produces  a 
conservative  estimate  (i.c.,  flows  that  are  larger  than  those  observed  during  historical  failures 
of  similar  blockages)  of  the  possible  peak  discharge  that  could  result  during  a  breaching  of 
the  Castle  Lake  Blockage. 

6)  The  effect  of  the  SRS  is  significant  in  reducing  the  peak  discharge  from  the  hypothetical 
failure  of  Castle  Lake.  The  initial  elevation  of  the  lake  prior  to  failure  also  affects  the 
magnitude  of  the  resulting  dambreak  discharge.  For  the  failure  and  bulking  scenario 
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recommended  by  the  Corps  of  Engineers,  the  SRS  reduces  the  peak  discharge  into  the  North 
Fork  Toutlc  River  by  85  percent  (from  695, 000  to  105,200  cfs)  for  full  lake  conditions.  If 
Castle  Lake  is  lowered  30  feet  prior  to  its  failure,  the  SRS  reduces  the  peak  flow  by  82 
percent  (from  352,500  to  62,(XK)  cfs).  If  Castle  Lake  is  lowered  by  60  feet,  the  SRS  reduces 
the  peak  flow  by  95  percent  (from  131,800  to  6,000  cfs). 

7)  The  amount  of  storage  the  SRS  can  provide  depends  on  how  full  of  sediment  it  is  when  a 
flood  event  occurs.  Under  the  present  "existing  conditions"  in  the  SRS  and  the  Cowlitz 
River,  all  of  the  Corps  recommended  flooding  scenarios  would  be  fully  contained  at  Castle 
Rock  and  Kelso  -  Longview.  The  resulting  flows  would  be  similar  to  a  100  year  flood  event 
in  the  Cowlitz  River.  If  the  SRS  were  full  of  sediment,  all  of  the  Corps  recommended 
flooding  scenarios  would  be  fully  contained  within  the  channel  at  Kelso  -  Longview,  but  not 
at  Castle  Rock  for  cither  lake  full  or  lake  lowered  30  feet  conditions. 

8)  None  of  the  hypothesized  breaching  and  bulking  scenarios  will  exceed  or  overtop  the  SRS  for 
either  "existing  conditions"  or  "full  conditions." 

9)  The  SRS  protects  communities  and  those  river  sections  downstream  from  it,  but  docs  not 
affect  the  areas  upstream  from  the  SRS.  Therefore,  there  is  no  protection  above  the  SRS 
from  a  hypothetical  failure  of  Castle  Lake.  Because  of  this,  and  the  possibility  that  a  failure 
of  Castle  Lake  would  greatly  reduce  the  sediment  storage  capacity  and  effective  life  of  the 
SRS,  additional  field  investigations  pertaining  to  the  geotechnical  stability  of  the  Castle  Lake 
blockage  and  continuous  monitoring  are  recommended. 


28 


5. 


References 


Costa,  John  E.  (1984)  "Physical  Geomorphology  of  Debris  Flows,"  in  Developments  and  Applications  of 
Geomorpholoev.  by  Costa,  J.E.  and  Fleishcr,  P.J.,  eds.,  Berlin,  Springer- Verlag,  pp.  269-317. 

Society  of  America,  Inc.,  Boulder,  CO. 

Costa,  John  E.  and  Schuster,  Robert  L.  (1986).  "A  Perspective  on  Landslide  Dams,"  in  Landslide  Dams: 
Processes,  Risk,  and  Mitigation,  by  Schuster,  Robert  L.,  American  Society  of  Civil  Engineers,  New  York. 
New  York. 

Costa,  John  E.  and  Schuster,  Robert  L.  (1988).  "The  Formation  and  Failure  of  Natural  Dams,"  in 
Geological  Society  of  America  Bulletin.  Vol.  100,  pp.  1054-1068,  July  1988. 

Fairchild,  Lee  H.  (1987).  "The  Importance  of  Lahar  Initiation  Processes,"  in  Debris  Flows/Avalanches: 
Process.  Recognition,  and  Mitigation,  by  Costa,  John  E.  and  Wieczorek,  Gerald  F.,  eds.,  The  Geological 

Frcad,  D.L.  (1989)  "BREACH:  An  Erosion  Model  for  Earthen  Dam  Failures,"  Users  Manual,  National 
Weather  Service,  Hydrologic  Research  Laboratory,  Silver  Spring,  Maryland. 

Frcad,  D.L.  (1989)  "DAMBRK:  The  NWS  Dam-Break  Flood  Forecasting  Model,"  Users  Manual, 

National  Weather  Service,  Office  of  Hydrology,  Silver  Spring,  Maryland. 

Froelich,  D.C.  (1987).  "Embankment-Dam  Breach  Parameters,"  in  Proceedings  of  the  1987  National 
Conference  on  Hydraulic  Engineering  at  Williamsburg,  VA,  American  Society  of  Civil  Engineers.  New 
York,  New  York. 

Laenen,  Antonius  and  Orzol,  L.  L.  (1987).  "Flood  Hazards  Along  the  Toutle  and  Cowlitz.  Rivers, 
Washington,  from  a  Hypothetical  Failure  of  Castle  Lake  Blockage,"  Water-Resources  Investigations 
Report  87-4055,  U.S.  Geological  Survey,  Portland,  OR. 

MacArthur,  Robert  C.  and  Hamilton,  Douglas  L.  (1989).  "Numerical  Simulations  of  Mudflows  Induced  by 
Hypothetical  Failure  of  the  Castle  Lake  Blockage  near  Mount  St.  Helens,  Washington,"  Projects  Report 
No.  89-04,  Hydrological  Engineering  Center,  Davis,  CA. 

MacDonald,  T.C.  and  Langridge-Monopolis,  J.  (1984).  "Breaching  Characteristics  of  Dam  Failures,"  in 
Journal  of  Hydraulic  Engineering.  American  Society  of  Civil  Engineers,  Vol.  110  (5),  pp.  567-586. 

Meyer,  D.F.  and  Dodge,  J.E.  (1988).  "Post-Eruption  Changes  in  Channel  Geometry  of  Streams  in  the 
Toutle  River  Drainage  Basin,  1983-85,  Mount  St.  Helens,  Washington,"  Open-File  Report  87-549,  U.S. 
Geological  Survey,  Vancouver,  WA. 

Schaefer,  M.G.  (1990).  Personal  communication  and  hand  written  notes.  Dam  Safety  Section,  Water 
Recourccs  Program,  Department  of  Ecology,  Olympia,  WA  98504-8711. 

U.S.  Army  Corps  of  Engineers,  Portland  District  (1990),  "Unpublished  Data:  Soil  Material  Properties  for 
the  Area  Below  Castle  Lake  to  the  N-l  Structure,"  collected  by  the  Geotechnical  Branch  of  the  Portland 
District,  July  1990,  Portland  OR. 

U.S.  Army  Corps  of  Engineers,  Portland  District  (1988).  "Castle  Lake:  Engineering  Analysis  and 
Alternative  Evaluation,"  Portland  OR. 

U.S.  Army  Corps  of  Engineers,  Portland  District  (1984).  "Mt.  St.  Helens  Cowlitz  and  Toutle  Rivers 
Sedimentation  Study,"  Portland,  OR. 


29 


U.S.  Army  Corps  of  Engineers,  Hydrologic  Engineering  Center  (1985).  "Toutle  River  Mudflow 
Investigation,"  Special  Projects  Report  No.  85-3,  Davis,  CA. 

U.S.  Department  of  the  Interior,  Bureau  of  Reclamation  (1973).  Design  of  Small  Danis,  a  Water 
Resources  Technical  Publication,  Washington,  D.C. 


30 


Appendices 


Appendix  A 

Annotated  Bibliography 

Prepared  by  Water  Engineering  and  Technology,  Inc., 

Fort  Collins,  CO 


X 

§ 


o 

*  -2  i 


3 

cr 

O 


D  -Q  D 
C  V 

m  is  D  o) 

3  a:  C 

0  4 1  — 

— -  — *  -* 


>-  "O  v» 

g  £  '3 


O  cfl 

3  <0 

CD  <D 


T3 

4) 


X) 

0) 

TJ 


~  T  I  £ 


«  -  £ 


^  CT> 

—  C  t- 

C  CD  O 

—  3=  > 

4)  °  ° 


(D  <D 

vi  -»  X) 

(TJ 

£  ?  * 

-  ■  E 


u> 


<0 


-O 
O  Q> 
«-  O 


4i 

*r  >* 

Q-  v> 


CD  —  _ 
CT>  Q>  O 

5  S  » 

M  c  £ 

jC  O  l_ 

0—0 

i  ts  I 


■S  c 
*-*  — 

S  (/) 


Hi  u 

5—3 

V)  T)  ■0 
O  <TJ  C 


CA 

O  (/) 

l_  (TJ 

Q.  flJ 

L. 

O  O 

(-  C 

(TJ  — 


V  3 
</>  — » 
CD  — 


& 


-  S 


0  —  0 


3  & 


I  i 


c  « 

4>  r-i  . 

O 

o  .  co 

</>  >  in 


-*  3  CA 

(TJ  V  »  • 

G>  OJ  OC  —  K) 

t~  C  —  nO 

CO  —  u  o  •— 

•  *-•  o  C  • 

E3*-3- 

(T)  O  <o  o  ro 

O  Of  3  U  «“ 


c/>  >v 
o  -  o 

§  >-  o 

3  —  v  O 

®  </>  <D  C 

TJ  >  U  C7> 

>*  —  o  c 

3:  o  to  uj 


V  ft 


* 

X  o 
o»  m 


>o 

N- 

S- 


Hi 

1 


•  ? 

>  <0 


oo—o 

—  3  U 

v*  O  O  3 

(TJ  —  —  — 


E 

<o 

o 


O  a 

L.  — 

-Q 


g  2 


?  a 

r  i 

C  CD 

o>  — 

i  5 


5 

*-> 

3 


?  . 

<TJ  U 


CJ  (A 
3 

-  i_ 


O 

QC 

s 


3 

—  CD 

—  c 


I  s 


o 

00 

o 


3 

—  o> 

—  c 

CD  — 

s  i 

O  X 


<0  "g 


1C 


a 

3 


A-2 


from  Hazard  Lake,  Yukon  Glaciology, 
Territory,  and  the  po.  3-21. 

proelem  of  f I ood 
magnitude  prediction. 


Rheology,  1988  Costa,  J.E.  Rheologic,  geomorphic  and  Flood  Presents  information  on  water  floods,  which  are  Newronian  fluids 

Myperconcentrated  sedimentologic  Geomorphology:  and  have  viscosities  that  are  unique  to  a  particula-  fluid 


c 

o 


<L»  -p 
p  6 
C  *-»  (A 

O  to  3 


>•  —  *  C 

O)  <u  .  o 

o  .xr  w>  (/>  • 

— '  O  TJ  Q. 

O  O  &  -Q  Kl 

-C  *  CO 

Q.  •  ID 

U'S  ,  - 

E  o  p  U  CQ 


o 

</> 

*—  o  </> 

8.0)0 
<0  O 


(/)  •—  O 

C  TJ  • 

0)  •  c 

*->  o  <u 

x  c  a 


C  3 

w  o 

C  -  t/i 

O 

•—  ro  a> 


i  °  §■ 

QUO) 


O  3  TJ  *-> 
—  -X  C  </) 


C 

o  — 

—  0) 

</)  c 

$  i  I  ? 

«o  *-  o  — 

<1>  13  — 

*-»  •  c  o 

s  ?  n 

o  o  o  i: 


floodplains. 


A- U 


St.  Louis, 
Missouri,  October 
26-30,  1981. 


X 

5 


O) 
L.  O 
Cj  — ' 

-C  o 

♦-»  L. 

«  p 

4)  > 

3  3: 


4) 

2  ? 


*-  a*  -n 

Qi  O  >  C 


g  £ 

•-  > 

4» 


i-  m  at 

IB  «-  C 

4>  O  •— 
w  -O  t- 

4»  <0  Q. 


w  a  j  w 


Cfl 


T>  41 
Q  L. 

6  3 

C  ^ 
O  <0 
•  —  H- 

y> 

2  § 

fti  T3 


ac  <o 

o  4> 

< 


nr 

11 


a>  — « 

x:  **-  o> 

»-  c 
^  ■— 
I  to  *-* 

•  4)  (A 

11  <-  <0 

OC  -fl  u 


* 

4> 


r  ? 

«  ■— 

i  i 


3 

? 

<0  — 

1 1 


? 


a  —  a  m- 


->  «*  C7> 

<a  <u  o 

c  o  —■*  o> 

o  —  o  c 

->«-•*“ 
w  c.  T>  u 

<0  a>  >  a 

JZ  </)  3C  00 


°  ®  s 

c/>  C 

O)  o  4* 

C  —  O  O 

C  •*- 
■n  n)  Qf  — 

y  i  l  3 

a  4i  « 

O  N.  *4-  U 

o  oo 

Jh  O  - 

o-  ■-  U 


E 

<0 

TJ 


• 

ii; 

4>  «  c 

-  w  5» 


*  o 
>*  2 
o> 

O  t_  sQ 

-»  a?  a 
o  a 

U  <0  ■ 

tj  a.  r\j 

o 

i  TJ  ^ 
4> 

15a 

> 

O  4» 

</>  C/>  IA 


>  • 
o>  o 
o  o 


c  to 

5  o  8 


C  t_ 

s  s 

c  2 


u 

Q 


O  -* 

-  > 
C  M 
•-  >- 
c-  4> 
>  oc 

N 

<0  "g 


00 

O 

o 


is 

4J 


*->  CO  ■— 

-  3  -  « 
t  .it 

<0  Q»  rtj  4> 
41  C  U 

m  JZ  §*  « 

■  ni 


A- 5 


4» 

g 


3 

<TJ 


Swiss  AIds. 


Failure  1982  Hagen,  V.K.  Re-evaluation  of  design  Paper  presented  at 

floods  and  dam  safety.  14th  International 


t_ 

o 


C  'P  N 

-  I  z 

"8  4J  J 


C  C/> 

O  c/> 

0 

c  i-  o 
0  0)1- 

tfl  o  5/ 
woe 


0  c 

C  —  JO  0> 

0  3  >o 

O.  «-*  O.  (M 

— •  <  5  -  C> 

0  O  <✓>  ID 

U  f-  I  X  W 


i-  0)  0 

-o  a  -q  <u  oo 

>  M  C  U  to 

:c  <  0  <  1- 


o  o 

0 

O  >.  • 

0  -3 

•  >  • 

(/>  «-  o 

•  3  ID 

D  GO  00 


XJ  •*- 

O  <0 
—  >- 
©  C/J  0 

3  <3  0 

*->  o  e 

a:  h-  3: 


©  i  « 

0  V) 
•4-^0 
O  — ' 

-  c 

>  0) 

0  0 

— 

0  S  0 

>  C/>  Q- 


0  w  — »  •  CO 

-C  O  O  *-»  O 

O  O  H-  X  «- 


->  c 
*  0 
-  »  y 


3  8  ?  3. 

3C  — I  0  (/) 


(X  V) 

c 

0  0 

1*1 

O  -  S 

U  GO  GO 


A-6 


Sabo  Works,  Japan  1976  Ikeya,  H.  Introduction  to  Sabo  The  Japanese  Sabo  Contains  an  introduction  to  tne  geology  and  natural  disasters  of 

Works-The  preservation  of  Association,  Japan.  Includes  history  of  Sabo  Works.  Describes  the  use  of 
land  against  sediment  Sankai-dc,  Co.,  Sabo  Works  for  sediment  control  in  channels  and  on  hillsides, 
disaster.  Tokyo,  168p.  T’-st  printed  in  English  in  'i>79. 


3 

O 


3 

(/> 


g 


-  P 

O  O 

>V  *” 

5  g 

L  <0 

3 

<0  o 
o 

<0  O 

co 

c  c 
•—  a* 
a  a> 

«-*  3 

c  *-< 

O  - 


Ji 


g 

g 


_  co 

O)  CO 

■O 


CO 

CO 

CO 

o 

CO 

c 

*-> 

<u 

o 

o> 

u 

a> 

o 

o 

CO 

(0 

l_ 

o 

3 

o 

a. 

00 

-C 

0) 

o 

o 

> 

CO 

a> 

g 

(0 

CO 

CO 

> 

L- 

(0 

<0 

3 

Q 

3C 

13 

CO 

6 

"O 

g 


c 

TJ 


1 1  t 

S(/>  O 
t/> 


P 

O' 


& 


«-  *o 

K\ 
C-  CO 


CO 

u 

flj 


o  o 

L. 

CO  <0 


&  i 


E  *- 
ft)  •— 
"O  — 

3 

to  O 
—  u 

g  g 

"O  <0 


*»  S#  "2 

O)  T3  <0 

C  ?  ' 

C  x:  oc 


*  r  ii 


jr  at 

?  *  £ 


S  8. 


a)  u  d 
—  <d  5. 


—  *  o  _ 

CO  CT>  O  C  - 

O  C  CO  -  OJ 

—  o» 

cue  c  •  • 

.c  a>  <o  uj  o  *- 

O  fc>  O  t  S- 

o  c  —  —  «- 

4-»  —  I—  —  •*  • 

o  o>  $  >0*0 

0)  C  g  o  in 

O  UJ  C  c_)  *—  •— 


(0 


0) 

>  .5 

3  E 

c  *«- 


-  g 

c  <u 

flj  </> 

i  g 


«  ° 
a  «  « 


U  ■  .£ 


C  -C  *♦- 

1  r  ° 
—  0(0 
i_  **-  <u 

8.  ~  *3 

2  3.  - 


CO 


(0 


—  (0  0> 

1_  ■—  (0 

a>  (o  >> 

(0  >*  — * 

-*  <0 

■  2  S 


g 


<0  co  m 

C  >  00 
—  Or 

«J  .V  *- 


M  (0 


o  «-  a 

—  Ol  Q. 

— '  *•  o  <u 

-  -  o  * 

CO 

c 


c 

-  _  —  <M 

o  -  - 

— «  £  O  ro  UJ 

ro  u  a>  o 

C  «J  C  —  — 


O  o* 

C  CO 


D  O  CT»  O  >  o  CM 

O  0»  C  E  —  O  K. 

“>  O  UJ  •<(_>*-*- 


c 

-  6 
*-»  (0 
m  Q 

L-  1~ 
O  OJ 
3 

c  o 


E 

to 

o 


o  — 
to 
c 


■—  a> 


CO  3 

5  =  g 

C  CO  CO 

—  H-  T3 


-  <u 
c  a> 
(0  (/> 


§  g 


CO 

o 


CO 

-  o 
c  — 
o  — 

CO  ~ 

c 

g  g 

“>  CT3 


>o 

N. 

o 


1 


CO  3 

a  — « 

«-* 

ro  rtJ 
U  U. 


.  &  * 

«-  Cl 
■  -O 

«-*  >  o 

t:  a  x 


A-7 


Construction,  Dec. 
1976,  pp.  C3-4S. 


CJ 

oc 

3 


i  - 

1  i 

CA  — • 

*4- 

11 

•2  g 

■— 

« 

Q-  (A 
I)  3  ^ 

CD  l-  t_ 

a»  <u  o 

>  o  3 

— *  ^ 

2  4»  O 

*—  CA  fc. 

D>  C  *■* 

rig 

O  (A  u 

ft) 

c  t_  c 

°  c  5 


£  -c  «) 

Q.  w 
3  •  *. 

l_  4~*  O 

ft)  (0  _D 


O  3  p» 
<D  ft)  C 
'♦-CO) 


*  6  a 

u  *-*  £ 

o  t a 

(A  o  n 

<U  u  C 

o  u  a> 


o  — 

ID  S 
ft) 


O  vt  ^ 
•-*-(/> 
o 

&  d  * 


a;  o 

-c  •«- 


o» 


00 
*♦-  o 

O  X- 


X  *-*  m 

-&  *  J 


-  «3  ^  _ 

L  OJ  « 

o>  >  o 

t-  0  —  u 

O  ft)  C  — 

'*-0  3  3: 


k_  u 

Is 


^  K- 
>  Q-  r>x> 

c  ?  o 


rO 

oo 

O 


-c  .  I\j 

jH 

ft)  *- 
o 

u.  • 
N  Z  & 

■o  i 

Q.  ^ 


o  g 

?  ? 
•  ID 
C 

o  c 
•—  o 


2  i  2 


V  (A 


*-*  o 

s  - 


c 

V  3 

--'  o 

CA  O 
3 


2  2 
~  ”2 


>» 

c 

c 

s*_ 

o 

CA 

c 

X*- 

u 

ft) 

o 

ft) 

CD 

o 

Oi 

c 

x: 

H— 

L. 

ID 

o 

p 

c 

X 

3 

u 

i 

V 

o 

CT 

L. 

ID 

p 

£ 

o 

L_ 

0) 

U 

o 

p 

4-> 

o 

*-* 

o 

CA 

i_ 

•4- 

‘4- 

'4- 

ft) 

■O 

l_ 

co 

ft) 

2 

fl) 

*D 

<D 

o 

> 

o 

(A 

<A 

2  2 


.2  g 

e  — 


%  s 

a  i 


S.  ; 

«  CD 

S  2 


f  I  o  ? 


o 

?  % 
1 0  T3 


^  O 
•—  CA 

3  U 


id  ft;  o» 

>  3  > 

<D  u  CO 


-Q 

ft> 

TJ 


— >  ID  O) 


C 

cu  ■— 

’  § 
TJ  T> 

<U  (A 

£  ii 


«D  C 
o>  - 

s  e 

i_  cu 


ft)  «-»  ►- 


-  £k 
l.  ft; 

a  -o 


<A 
a> 

CA  <A 
3  ft) 

■o 


■Z  2 


& 


O  3 
ID  (A 
as  3 

6 

i—  c 
01 

V  * 

(a  m 

ID 

(A  • 

—  > 


c  c 

””  5 


S  21 

^  O 

^  *CA 


(A  *-> 
(A  {A 
g  CD 


S  a 


K. 

O' 


»0  <A 
>•  -C  ^ 
-Q  - 


2 


«j  o> 

£  QJ 
>  O 

i  2 


o  — 
o  o 

ID  > 


11 

*  o 


ID 

3  CD 

§>  i 
i  2 


C 

a 


3  C 

CA  O 
3  t- 

•  > 
m  C 
X  Hi 


o 

ft» 

o 


A- 8 


trigger  a  downvalley  mass  movement  were  small  (8-15  rim  for  one-four 
intensity). 


fO 


c 

I 

w 

l_ 

4) 

5 


A- 9 


Landslide  Dam  1960  Knight,  D.K.  Stability  of  slide  dam  Report  on  Flood 

and  Bennett,  and  recommendations  on  Emergency,  Madison 


£ 


4>  — ' 

>  r- 

43  CL 

*0  c/3 


o  » 
~o 

2  o 


43  43 

x:  •*- 

**  "S  • 

O  o 


«  8 
•—  <0 
t- 

n  c 

43  O 
Q  ■" 
*-» 

?  2 
flj  C 


C  cn 
■—  c 

?  c 


c  o  m 

3  C 

O  crt  to 

U  *-<  <J 

u  — »  u» 

03  3 

C/3  O) 

-g  43  c 


C/3  — • 

C  .8 
43  43 

1/3  TJ 
43 


* 


o 

8.  iff 


to  c 

flO 

<D  w 


Y  = 


•—  4-* 

*-•  3 

O 


to 

4-* 

a 


2  S  .? 
g  °  2 
°  ■*-  2 

^  C  03 
4)  C 

§■  1 
CO  — * 

-S  « 

13)  — 

?  „* 


<U  XI 
TJ  4; 
2  O 


43  4-«  *-» 

4-*  03 

1  -  O 

-  3 

o  -a 

if  JC  c 

<U  c/3  CD 


(0 

c 

o 

(D  0) 

O  03 
—  0) 

03  H- 
O  O  X 


43 

CD 

• 

to 

to 

TJ 

to 

a 

c 

c_ 

W 

43 

O 

03 

> 

-* 

c 

(D 

LU 

OH 

a 

3) 

2L 


4)  L 

JC  C3 

O  03 

c  -  c 

03  O  43 

^  CO  w 

~  4» 


II2 


C  43 
m  o 
c 

c  <0 
°  u 

c/3  CO  4-* 

<n  «)  4-  to 

EC4IC- 
U.  03  X  .a  X 

o  —  <0  4/  <D  *■* 

u.  (/)  -1  TJ  I  I 


2  ? 
3  «J 

■8  o 

*  a  ? 
to  to  co 

4->  to  w- 
m  X. 

X  if* 
if  w 

X  03  § 

T)  t>  t- 
03  •— 
to 

*-<  x  t/3 

C/i  — •  *-» 

c  <c  c 

3  C  4l 


if 

to 

c  « 

—  43 

a.  *-* 

if 


TJ  4/ 

a>  o 

TJ  O 

>  C 
o  o 

a  ~ 


I  i 


£ 


43  -C 
§  .2 
<0  ”§ 


w 

$1 


c/»  <0 

C  -C  — 

—  o 

CD  •  to 

♦-*  C  43 

C  43  u 

O  fi.  H 

O  O  *-• 


—  C/1 


X  -O 
in  * 
to 

?  § 

CD 

03 

*—  C 
4/  ■— 
>  *-» 
•—  CO 
OC 

X 

4l  43 
01 

3  t 


fz  to 

r  ? 

o  _ 
z  o 
c 

c  o 
o  o 


*-  43 

c 

C/3  C  X 

ro  *-■ 
i _  .c  c 

X>  O  ~ 

Qj  CD 

■Q  4»  *-» 

-C  C- 
4-/  to 


s  ~ 
fe.  ? 


3  g 
y  5 


"43  13 

TJ  c 


o 

43 

43 

C_ 

9 

TJ 


E  •- 
CD 

43  3 


O 

43 

O  X 

•  5 

co  e 


—  -  (0 
M 

CD  05 
— «  (_ 

3  ? 


2  5  8. 

<0  C  4l 

3  —  Of 

X 
*  o 

C  — *  — ' 

O  CD  *4-. 

4^  c  6 

O)  O  CO 

C  —  43 


TJ 

43 

C 

o 


X 

T3 

CD 

43 


■3 


43  O  ? 
C33  n  o 

u  — ■ 

It  TJ  4- 
JZ  — 
to  3  43 


o 

XI 


8  CD 
X  *4- 


O 

03 

o 

o 

4> 


X 

43 

> 


CO  — 

3 

if  O 
TJ  cj 

s  ? 

CD  CD 
^  43 

*8  - 


S  8 

C  0) 
43  C 
CD  CO 


C 

<U  — 

c  o 
43  M 


o 

CO 

O 


->  43 

C 

5  43  X- 


C  *-  o  to 

X  .  g  ^  C  ^ 

5  2  i,  8  t  g 

j  I  O  «»  «/>  W 


A- 1 0 


River,  'Irod  Hazard  a  hypothetical  failure  of  Investigations  is  assumed  to  increase  by  a  factor  of  5  from  Castle  Lare  to  N-1  dam. 

Castle  Lake  blockage.  Report  87-105j.  The  hypothetical  hyperconcentrated  flow  is  routed  downstream, 

suoerimposed  on  normal  winter  flood  flows  by  use  of  a  one-dimensional, 
unsteady 'St ate  numerical  streamtlow  simulation  node l . 
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Dam  Break,  Dam  1980  Lana,  L.F.  Evaluation  of  selected  Geological  Survey 

Failure,  Modeling  aam-break  floodwave  Uater-Resojrces 

models  by  using  field  Investigation,  v. 


Mt.  St.  Helens,  1983  lisle,  T.E.,  Stream  channel  Proceedings  of  the  Describes  wide  variations  in  stream  channel  response  to  the  1980 

Geomorphology,  lehre,  A.K.,  adjustments  after  the  Symposium  on  eruption  of  Mt.  St.  Helens.  Cross  sections  and  longitudinal 

Channel  Ma-tinson,  1980  Mt.  St.  Helens  Erosion  Control  in  profiles  in  all  major  channels  that  drain  the  blast  area  were 
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Mudflow,  Toutle  1985  MacArthur,  Tootle  River  mudflow  Special  Projects  Presents  the  results  of  a  study  whose  primary  purpose  was  to 

River,  Sediment  R.C.,  Hamilton,  investigation.  Report  No.  85-3,  develop  and  apply  a  one- dimensions l  unsteady  mudflow  routing 

Retention  Structure  O.L.,  and  prepared  for  U.S.  model  to  route  estimated  mudflow  events  through  the  reservoir  and  its 
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Mt.  St.  Helens,  1981  Meier,  M.F.,  Hydrologic  effects  of  Mt.  EOS,  Transactions 

Volcanic  Dam,  Carpenter,  St.  Helens'  1980  of  the  American 

Overtopping,  Dam  P.J.,  ana  eruption.  Geophysical  Union 
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Oebris  Flow  1988  Okunishi,  K.,  Hydrological  controls  of  Hydrolog;cal  Discusses  ♦our  years  of  hydrological  observations  as  related  to 

Suwa,  H.,  and  erosion  and  sediment  Sciences  Journal,  the  triggering  of  debris  flows  in  the  Kami -kamihori  Valley  in  the 
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NUMERICAL  SIMULATION  OF  MUDFLOWS  FROM  THE  HYPOTHETICAL  FAILURE  OF 
A  DEBRIS  BLOCKAGE  LAKE  BELOW  MOUNT  ST.  HELENS,  WA 

Robert  C  MacArthur1,  M-ASCE,  Douglas  L  Hamilton2,  M-ASCE,  and 
Ronald  C  Mason1,  MASCE 


Abstract 

This  paper  evaluates  the  characteristics  of  mudflow  events  resulting  from  the  hypothetical 
breaching  of  a  debris  blockage  dam  using  a  variety  of  lake  levels  and  impounded  water  volumes 
for  the  initial  breach  conditions.  A  one-dimensional  (Petrov-  Galerkin  finite  element)  unsteady 
mudflow  routing  model  is  used  to  simulate  the  movement  of  the  dam  break-induced  mudflow 
events  down  valley  through  the  Corps  of  Engineers*  Sediment  Retention  Structure 

Introduction 

The  May  18,  1980  eruption  of  Mount  SL  Helens,  WA,  produced  a  debris  avalanche  which 
t  -»wed  down  the  North  Fork  Toutle  River  damming  several  tribut  .y  streams.  The  blockage  at 
the  confluence  of  South  Fork  Castle  Creek  and  Castle  Greek  produced  a  natural  ibris  dam 
approximately  190  high.  Figure  1  shows  the  general  study  area  near  Mount  St.  Helens  and  the 
location  of  Castle  Lake.  Snow  melt  and  runoff  waters  captured  behind  the  blockage  quickly 
formed  a  lake.  To  prevent  overtopping  and  a  potentially  catastrophic  failure  of  the  blockage 
retaining  Castle  Lake,  the  U.S.  Army  Corps  of  Engineers  (COE)  constructed  an  SPF  spillway  at 
the  eastern  end  of  the  blockage  to  stabilize  the  lake  at  elevation  2,577  feet  MSL  Recent  studies 
by  the  U.S.  Geological  Survey  (USGS)  indicate  that  'the  blockage  is  potentially  unstable  against 
failure  from  piping  due  to  heave  and  internal  erosion  when  groundwater  levels  are  seasonally  high' 
and  that  an  earthquake  of  6.8  or  greater  might  initiate  such  a  failure  (Laenen  and  Orzol,  1987). 
If  the  Castle  Lake  blockage  were  to  fail  rapidly  by  the  mechanism  suggested  by  the  USGS, 
approximately  18,500  acre-feet  (AF)  of  stored  water  in  the  lake  could  create  a  mudflow  flood  event 
in  the  North  Fork  Toutle  River.  The  USGS  (Laenen  and  Orzol,  1987)  estimates  that  an  event 
of  this  nature  could  result  in  a  peak  discharge  of  2,100,000  cfs  at  the  Corps’  N-l  debris  retention 
dam  twelve  miles  downstream  from  Castle  Lake  (see  Figure  1). 

In  May  1982,  President  Reagan  directed  the  Corps  of  Engineers  to  prepare  a  comprehensive 


Hydraulic  Engineer,  U.S.  Army  Corps  of  Engineers,  Hydrologic  Engineering  Center,  609  2nd  St,  Davis,  CA  95616,  U.SA. 

2  Principal,  RIVERTECH,  Inc,  23332  Mill  Creek  Drive,  Laguna  Hills,  CA  92653,  U.S.A. 

3  Chief,  Hydrologic  Section,  U.S.  Army  Corps  of  Engineers,  Portland  District,  P.O.  Bon  2946,  Portland,  OR  97208  U.SA. 
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plan  for  long-term  flood  control  and  navigation  maintenance  in  the  wake  of  the  Mount  St.  Helens 
eruption-  A  major  component  of  the  resulting  plan  is  the  recently  constructed  Sediment  Retention 
Structure  (SRS).  The  primary  function  of  the  $56.5  million  dam  is  to  trap  the  huge  amounts  of 
sediment  expected  to  continue  to  move  down  the  North  Fork  Toutle  River.  The  CHS  was 
designed  to  capture  runoff-induced  sediment  from  the  blast  zone,  thus  preventing  sediment 
deposition  and  reduced  flood  routing  and  navigation  capacities  in  the  Cowlitz  and  Columbia 
Rivers.  Failure  ot  the  Castle  Lake  blockage  resulting  in  the  possible  occurrence  of  a  mudflow 
event  could  jeopardize  the  safety  and  performance  of  the  SRS. 

The  purpose  of  this  study  is  to  evaluate  the  hydraulic  characteristics  of  mudflow  events 
resulting  from  the  hypothetical  failure  of  Castle  Lake  and  to  examine  the  ability  of  the  SRS  to 
capture  and  pass  such  events  through  its  spillway  for  various  initial  conditions  at  Castle  Lake  and 
in  the  SRS.  More  spedflcally,  the  study  is  to:  (1)  determine  if  flows  will  exceed  the  present 
spillway  capacity  at  the  SRS,  (2)  determine  if  the  SRS  will  be  overtopped,  (3)  estimate  how  the 
peak  discharge  downstream  from  the  SRS  will  be  affected  by  the  presence  of  the  SRS,  (4)  evaluate 
the  effects  on  the  resulting  routed  mudflow  hydrograph  due  to  lowering  the  initial  Castle  Lake 
levels  at  the  time  of  breaching,  and  (5)  evaluate  the  performance  of  the  SRS  during  these  various 
events  when  the  SRS  is  empty  with  (a)  'existing  conditions',  (b)  'half  full"  of  accumulated 
sediment  and  debris,  or  (c)  'completely  full*  of  sediment  deposits  up  to  the  spillway  crest 

Approach 

The  Hydrologic  Engineering  Center  (HEC)  applied  the  one-dimensional  (Petrov-  Galerkin 
finite  element)  unsteady  mudflow  routing  model  (MacArthur,  et  aL,  1988)  to  route  several 
hypothetical  mudflow  events  from  Castle  Lake  to  the  SRS.  The  Mudflow  Model  was  modified  to 
incorporate  the  400  foot  wide  spillway  at  the  SRS  as  well  as  the  possibility  of  overtopping  the 
structure.  Therefore,  the  compute  outflow  hydrographs  downstream  from  the  SRS  include  effects 
due  to  storage  (ponding)  inside  the  SRS,  flow  through  the  spillway,  and  overtopping  of  the  SRS. 
The  upstream  boundary  of  the  modeling  reach  was  established  at  the  N-l  structure  to  correspond 
with  the  location  of  the  USGS  developed  mudflow  hydrographs  for  various  initial  lake  elevations 
prior  to  breaching.  The  routing  reach  is  approximately  6.1  miles  with  an  average  bed  slope  of 
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0.0093ft/ft.  HEC  applied  the  same  dcbulking  mechanism  downstream  lrom  the  N-l  structure  that 
the  USGS  prescribed  (Lcancn,  et  al.,  1990). 

HEC  used  topographic  data  and  measured  cross-section  information  prepared  by  the  USGS 
(Laenen  and  Orzol,  1987)  and  by  the  Portland  District  Corps  of  Engineers  to 
depict  the  natural  valley  geometry  of  the  mudflow  routing  reach  from  the  N-l  structure  to  the 
SRS.  Fluid  properties  used  to  describe  the  rheological  characteristics  of  the  USGS’  hypothetical 
mudflow  were  obtained  from  Major  (1984).  He  reported  fluid  viscosities  from  the  1980  Mount 
St  Helens  eruption  ranging  from  6  to  100  lb.-sec/sq.ft,  yield  strengths  from  2  to  31  lb./sq.ft,  and 
unit  weights  from  100  to  125  lb./cu.ft  For  the  purposes  of  this  investigation,  HEC  used  the 
following  constant  fluid  properties:  Viscosity  =  6.0  lb.-  sec/sq.ft,  yield  strength  =  2.0  lb./sq.ft., 
and  unit  weight  =  110  lbVcu.ft  These  fluid  properties  were  chosen  to  give  the  most  conservative 
results  (e.g.,  the  greatest  velocities  and  the  least  attenuation).  Sensitivity  and  model  validation 
studies  previously  conducted  by  MacArthur,  et  al.,  (1985  and  1987)  support  this  reasoning. 

Three  different  valley  and  SRS  geometry  scenarios  were  evaluated.  First,  the  SRS  was  assumed 
to  be  empty  of  water  with  the  present  (existing)  dry  reservoir  bottom  at  elevation  870’  NGVD  and 
the  valley  and  channel  upstream  from  the  SRS  represented  by  the  Corps  of  Engineers 
photogrametric  cross-sections.  These  conditions  are  referred  to  as  ’existing  conditions*  throughout 
this  paper.  The  second  scenario  assumes  that  the  SRS  is  "half  full’  of  sediment  deposits  and  t..c 
channel  bed  in  the  North  Fork  Toutle  River  immediately  upstream  from  the  SRS  is  full  of 
sediment  with  a  bed  slope  of  0.006  ft/ft  (approximately  half  the  original  bed  slope).  The  third 
scenario  assumes  the  SRS  is  ’full"  of  sediment  deposits  up  to  the  spillway  crest  elevation  of  940’ 
NGVD  and  the  bed  slope  in  the  channel  upstream  from  the  reservoir  is  at  0.006  ft /ft 

In  order  to  evaluate  the  downstream  effects  of  lowering  the  initial  lake  levels  behind  the  Castle 
Lake  blockage,  five  different  initial  lake  surface  elevations  were  used  to  develop  five  different  dam 
break  hydrographs  (Laenen,  et  al.,  1990  Unpublished  report).  Those  five  dam  break  hydrographs 
were  bulked  up  to  reflect  the  likely  entrainment  of  sediment  and  debris  materials  (Laenen  and 
Orzol,  1987).  Each  of  the  five  hydrographs  was  routed  from  the  upstream  boundary  of  the 
modeling  reach  at  the  N-l  structure,  downstream  to  the  SRS  for  each  of  the  three  initial  SRS 
sedimentation  conditions  -  (a)  existing,  (b)  half  full  and  (c)  full.  Figures  2,  3  and  4  show  the  five 
bulked  mudflow  hydrographs  at  the  N-  1  resulting  from  each  of  the  five  starting  lake  elevations 
and  the  computed  outflow  hydrographs  from  the  SRS  for  each  of  the  three  initial  SRS  conditions. 

The  mudflow  model  was  modified  to  include  the  effects  of  the  400  foot  wide  spillway  at  the 
SRS  as  well  as  the  possibility  of  overtopping  the  structure.  The  modelers  assumed  that  the  six 
rows  of  3-foot  diameter  conduits  in  the  SRS  would  clog  and  become  inoperative  during  events  of 
these  magnitudes.  This  is  a  reasonable  and  conservative  assumption.  Therefore,  routed  mudflow 
hydrographs  downstream  from  the  SRS  include  effects  due  to  storage  (ponding)  inside  the  SRS, 
flow  through  the  spillway,  and  overtopping  of  the  SRS. 

Results  and  Discussion 

Computed  travel  times  for  the  leading  edge  of  the  mudflow  bores  are  summarized  in  Table  1. 
They  range  from  a  low  of  25  minutes  for  existing  SRS  conditions  with  2.10  Mcfs  inflow  to  a  high 
of  29  minutes  for  a  full  SRS  with  the  same  inflow.  Therefore,  mudflows  resulting  from  a  ’full 
lake*  breach  and  entering  an  ’existing’  SRS,  or  a  "full"  SRS  will  move  downvalley  from  the  N-l 
at  speeds  of  approximately  21.5  ft/sec  and  18.5  ft/sec,  respectively.  Lowering  the  initial  Castle  Lake 
elevations  by  60  feet  increases  the  computed  travel  times  t  51  minutes  and  65  minutes, 
respectively  (flow  velocities  of  10.5  and  83  ft/sec).  Table  1  and  t  ,ures  2  and  3  summarize  the 
five  different  inflow  and  routed  outflow  hydrographs  for  the  three  different  initial  SRS  conditions 
(existing,  half  full,  and  full),  respectively.  The  peak  discharges  at  the  upstream  model  boundary 
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Table  1  Castle  I^ikc  Dam  Breach  Routing  Results 


fl) 

(2) 

(3) 

(4) 

(5) 

(6) 

Initial 

Peak 

Peak 

Peak 

Average 

Castle 

Inflow 

Outflow 

Stage 

Travel 

Wave 

Run 

SRS 

Lake 

at  N-l 

at  SRS 

at  SRS 

Time 

Velocity 

No. 

ConditionEIcv(ft) 

(cfc) 

(cfc) 

(ft) 

(min) 

(ft/sec) 

1 

Existing 

2577 

2.10  M 

194,000 

974 

25 

21.5 

2 

Existing 

2562 

1.41  M 

71,000 

957 

27 

19.9 

3 

Existing 

2547 

1.18  M 

0 

938 

31 

173 

4 

Existing 

2532 

0.728  M 

0 

918 

39 

133 

5 

Existing 

2517 

0.437  M 

0 

902 

51 

103 

6 

Half 

2S77 

2.10  M 

244,000 

987 

24 

20.6 

7 

Half 

2542 

1.41  M 

189,000 

975 

->9 

183 

8 

Half 

2547 

1.18  M 

108,000 

942 

33 

143 

9 

Half 

2532 

0.728  M 

0 

920 

39 

133 

10 

Half 

2517 

0.437  M 

0 

908 

51 

103 

It 

Full 

2577 

2.10  M 

400,000 

1,009 

29 

183 

12 

Full 

2542 

1.41  M 

411,000 

1,003 

32 

163 

13 

Full 

2547 

1.18  M 

315,000 

995 

47 

11.4 

14 

Full 

2532 

0.728  M 

249,000 

985 

49 

11.0 

15 

Full 

2517 

0.437  M 

230,000 

982 

45 

83 

(1)  Starting  take  water  surface  elevations  assumed  by  the  USGS  (Laenen,  ct  at.,  1990) 

(2)  Peak  inflow  at  the  N-l  structure  from  the  USGS  (Laenen,  et  at  1990) 

(3)  Mudflow  routing  results  at  the  SRS,  this  study 

(4)  Computed  peak  stages,  this  study 

(5)  Estimated  flood  wave  travel  time,  this  study 

(6)  Average  tip  velocity  from  N-l  to  SRS,  this  study 


(at  the  N-l  structure)  range  from  2.10  Mels  for  Castle  Lake  full  at  elevation  2,577  ft  MSL  to 
0.437  Mcfc,  assuming  Castle  Lake  has  been  lowered  to  elevation  2^17  fL  MSL  prior  to  breaching. 
The  maximum  computed  outflow  from  the  SRS  using  full  Castle  Lake  starting  conditions  ranges 
from  196,000  cfc  to  600,000  cfc  for  'existing"  and  "full"  SRS  conditions.  Results  presented  in  Table 
1  and  Figures  2  and  3  show  that  mudflow  events  for  run  numbers  1  through  10,  13,  14,  and  15 
are  contained  within  the  SRS  and  do  not  overtop  the  structure  or  exceed  the  designed  spillway 
capacity  of  340,000  cfc.  Mudflow  events  exceed  the  capacity  of  the  spillway  and  overtop  the  dam 
for  runs  11  and  12  only.  The  peak  outflows  from  the  SRS  under  these  conditions  are  600,000  cfc 
and  411,000  cfc  respectively. 

Therefore,  under  "existing  conditions’  and  the  initial  lake  level  at  2577,  the  SRS  reduces  the 
peak  discharge  in  the  North  Fork  Toutle  River  by  90%  percent  (from  2.10  Mcfc  to  196,000  cfc). 
For  "half  full"  conditions  the  peak  outflow  from  the  SRS  is  reduced  by  87%  percent  and  by  71% 
percent  for  "full  conditions."  All  five  mudflow  hydrographs  for  different  initial  lake  levels  are  fully 
contained  within  the  SRS  and  Spillway  without  overtopping  for  "existing"  and  "half  full"  conditions. 
Overtopping  occurs  for  initial  lake  levels  above  2547T  MSL  and  "full*  SRS  conditions.  For  the 
worst  case  conditions  (SRS  full  and  initial  lake  elevation  at  2577’  MSL)  the  peak  outflow  is 
600,000  cfc.  The  SRS  obviously  reduces  the  peak  discharge  in  the  North  Fork  Toutle  River 
downstream  from  the  SRS  for  all  three  initial  SRS  infill  conditions.  Lowering  the  initial  Castle 
Lake  elevations  at  the  time  of  the  assumed  breach  also  reduces  peak  flows  entering  and  leaving 
the  SRS.  HEC  incorporates  into  all  these  results  the  estimated  volume  reduction  due  to 
debulking  of  the  flows  below  the  N-l  as  prescribed  by  Laenen  and  Oral  (1987). 
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Figure  2  Inflow  and  Outflow  Hydrographs  for  SRS  "Existing"  Condition. 


Figure  3  Inflow  and  Outflow  Hydrographs  for  SRS  "Full"  Condition 


Summary  and  Conclusions 

The  one-dimensional  (Petrov-Galerkin  finite  element)  unsteady  mudflow  routing  model 
was  used  to  route  the  five  different  hypothetical  mudflow  events  from  the  Corps’  N-l 
structure  to  the  SRS  for  three  beginning-of-event  geometry  scenarios  (SRS  initial 
conditions).  For  "existing  conditions"  in  the  SRS  the  maximum  (lake  full)  mudflow 
hydrograph  is  reduced  by  90  percent  and  the  SRS  is  not  overtopped.  If  the  SRS  is  "half 
full"  of  sediment  deposits  when  the  maximum  (lake  full)  mudflow  occurs,  the  peak 
mudflow  hydrograph  is  reduced  by  80  percent  and  the  SRS  is  not  overtopped.  However, 
the  peak  stage  comes  within  10  feet  of  the  top  of  the  dam.  If  the  SRS  is  initially  "full" 
of  sediment  to  the  crest  elevation  of  the  spillway,  the  dam  is  overtopped  for  mudflow 
events  where  the  assumed  initial  Castle  Lake  water  elevations  are  higher  than  2547  MSL.  The 
following  conclusions  are  made  based  on  the  results  presented  in  this  report: 

l.The  SRS  reduces  flows  in  the  downstream  reaches  of  the  North  Fork  Toutle  River,  even  for 
a  maximum  mudflow  event  of  the  magnitude  estimated  by  the  USGS  (peak  Q  =  2,100,000  cfs). 
Maximum  peak  mudflows  are  reduced  approximately  90%,  87%  and  71%  if  the  SRS  is  in  its 
existing,  half  full  and  full  condition,  respectively. 

2-For  existing  and  half  full  conditions,  no  overtopping  occurs  at  the  SRS  and  the  peak 
discharge  into  the  downstream  reach  is  reduced  to  196,000  cfs  and  266,000  els,  respectively. 

3. For  full  SRS  conditions,  the  dam  is  overtopped  for  those  mudflows  that  occur  with  initial 
Castle  Lake  elevations  greater  than  2547  MSL  The  maximum  depth  of  overtopping  is 
approximately  10  feet. 

4. Reduction  of  the  initial  Castle  Lake  levels  significantly  reduces  the  magnitude  of  the  resulting 
dam  breach-induced  mudflow. 

5. AdditionaI  economic  analyses  are  necessary  to  evaluate  the  cost/benefit  characteristics  of 
constructing  mitigative  measures  for  reducing  the  initial  lake  level.  Studies  are  being  conducted 
by  the  U.S.  Army  Corps  of  Engineers,  the  U.  S.  Forest  Service  and  the  U.S.  Geological  Survey 
to  better  determine  the  "most  likely  dam  breach  scenario"  and  the  most  effective  way  to  insure 
safe  operation  of  the  Sediment  Retention  Structure. 
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Appendix  C 

Routing  Results  for  the  Different  Breaching  and 
Bulking  Scenarios  that  were  Considered 


The  Following  Plots  are  Results  from  the  Analysis  of  the  Breaching  and  Bulking 
Scenario  Represented  by  a  Piping  Failure  with  an  Ultimate  Bulking  Factor  of  2.50. 
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The  Following  Plots  are  Results  from  the  Analysis  of  the 
Breaching  and  Bulking  Scenario  Represented  by  a  Piping  Failure  with  an  Ultimate 

Bulking  Factor  of  3.32. 

This  Scenario  is  Referred  to  as  the  HEC  or  Corps  of  Engineers  Recommended 
Scenario  for  Breaching  and  Bulking. 


0-  I  ! 


HEC  SCENARIO,  SRS  TO  COLUMBIA,  BF  =  3. 33,  LAKE  FULL,  SRS  FULL 


AT  MILE  33. SI 
AT  MILE  47. SI 


COLUMBIA.  BB=2.32,  LAKE  FULL,  EXISTING  CONDITIONS  AT  SRS 


BF=3.32 


POUTED  RESULTS  AT  MILE  6.0  f\ - ROUTED  RESULTS  AT  MILE  12.86 

ROUTED  RESULTS  AT  MILE  1.87  <0, - ROUTED  RESULTS  AT  MILE  15.71 

ROUTED  RESULTS  AT  MILE  5.31 

ROUTED  PESULTS  AT  MILE  10.19 


60',  BF= 


EC  SCENARIO,  SRS  TO  COLUMBIA,  BF=C.32,  LAKE  LOWERED  60',  SRS  FULL 


ROUTED  RESULTS  AT  MILE  33 
ROUTED  RESULTS  AT  MILE  47 


HEC  SCENARIO,  SRS  TO  COLUMBIA,  BF=3.32.  LAKE  LOWERED  60\  EXIST.  COND.  AT  SRS 


The  Following  Plots  are  Results  from  the  Analysis  of  the 
Breaching  and  Bulking  Scenario  Represented  by  a  "Heave  Type"  Failure  Using  an 
HEC  Bulking  Procedure  and  an  Ultimate  Bulking  Factor  of  4.46. 
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